Distributed Source Coding using Abelian Group Codes 



Dinesh Krithivasan and S. Sandeep Pradhan* 
Department of Electrical Engineering and Computer Science, 
University of Michigan, Ann Arbor, MI 48109, USA 

email: dineshkOumich . edu , pradhanvOeecs . umich . edu 

August 19, 2008 

Abstract 

In this work, we consider a distributed source coding problem with a joint distortion criterion depending 
on the sources and the reconstruction. This includes as a special case the problem of computing a function 
of the sources to within some distortion and also the classic Slepian-Wolf problem [62] , Berger-Tung problem 
[46], Wyner-Ziv problem [44], Yeung-Berger problem [47] and the Ahlswede-Korner-Wyner problem [42,63]. 
While the prevalent trend in information theory has been to prove achievability results using Shannon's 
random coding arguments, using structured random codes offer rate gains over unstructured random codes 
for many problems. Motivated by this, we present a new achievable rate-distortion region (an inner bound to 
the performance limit) for this problem for discrete memory less sources based on "good" structured random 
nested codes built over abelian groups. We demonstrate rate gains for this problem over traditional coding 
schemes using random unstructured codes. For certain sources and distortion functions, the new rate region 
is strictly bigger than the Berger-Tung rate region, which has been the best known achievable rate region 
for this problem till now. Further, there is no known unstructured random coding scheme that achieves 
these rate gains. Achievable performance limits for single-user source coding using abelian group codes are 
also obtained as parts of the proof of the main coding theorem. As a corollary, we also prove that nested 
linear codes achieve the Shannon rate-distortion bound in the single-user setting. Note that while group 
codes retain some structure, they are more general than linear codes which can only be built over finite fields 
which are known to exist only for certain sizes. 

1 Introduction 

A large number of problems in multi-user information theory fall under the general setup of distributed source 
coding. The most general framework for a distributed source coding problem consists of a set of encoders 
which observe different correlated components of a vector source and communicate their quantized observations 
to a central decoder through a rate-constrained noiseless communication link. The decoder is interested in 
reconstructing these observations or some function of them to within some distortion as measured by a fidelity 
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criterion. The goal is to obtain a computable single-letter characterization of the performance limits measured 
by the rates of transmission and the distortions achieved. Such a formulation finds wide applications in many 
areas of communications such as sensor networks and distributed computing. 

There is a vast body of work that addresses this problem, and solutions have been obtained in a variety of 
special cases [5-7, 11, 39, 42, 44, 45, 47, 50, 630. All of the existing works use the following canonical encoding 
strategy. Each encoder has two operations, implemented sequentially, each of which is a many-to-one mapping: 
(a) quantization and (b) binning. In quantization, typically, neighboring source sequences are assigned a code- 
word, whereas in binning, a widely separated set of codewords is assigned a single index which is transmitted 
over a noise-free channel to the decoder. The decoder looks for the most likely tuple of codewords, one from 
each source, and then obtains a reconstruction as a function of this tuple of codewords. 

In most of these works, existence of good encoders and decoder is shown by using random vector quantization 
followed by random independent binning of the quantizer codebooks. The best known inner bound to the 
performance limit that uses this approach is the Berger-Tung inner bound. It has been shown in the literature 
that this is optimal in several cases. The work of Korner and Marton [39], however, is an exception and 
looks at a special case of the problem involving a pair of doubly symmetric binary sources and near lossless 
reconstruction of the sample-wise logical XOR function of the source sequences. They considered an encoding 
strategy where the first operation is an identity transformation. For the second operation, they consider random 
structured binning of the spaces of source sequences and show optimality. Further, the binning of two spaces 
is done in a "correlated" fashion using a binary linear code. In the present paper, we build on this work, and 
present a new achievable rate region for the general distributed source coding problem and demonstrate an 
encoding scheme that achieves this rate region by using random coding on structured code ensembles. In this 
approach, we consider the case where the sources are stationary discrete memory less and the reconstruction is 
with respect to a general single-letter fidelity criterion. The novelty of our approach lies in an unified treatment 
of the problem that works for any arbitrary function that the decoder is interested in reconstructing. Further, 
our approach relies on the use of nested group codes for encoding. The binning operation of the encoders are 
done in a "correlated" manner as dictated by these structured codes. This use of "structured quantization 
followed by correlated binning" is in contrast to the more prevalent "quantization using random codes followed 
by independent binning" in distributed source coding. This approach unifies all the known results in distributed 
source coding such as the Slepian-Wolf problem [62], Korner-Marton problem [39], Wyner-Ahlswede-Korner 
problem [42,63], Wyner-Ziv problem [44], Yeung-Berger problem [47] and Berger-Tung problem [46], under a 
single framework while recovering their respective rate regions. Moreover, this approach performs strictly better 
than the standard Berger-Tung based approach for certain source distributions. As a corollary, we show that 
nested linear codes can achieve the Shannon rate-distortion function in the single source point-to-point setting. 
A similar correlated binning strategy for reconstructing linear functions of jointly Gaussian sources with mean 
squared error criterion was presented in [54]. The present work develops a similar framework based on group 
codes. 

This rate region is developed using the following two new ideas. First, we use the fact that any abelian 
group is isomorphic to the direct sum of primary cyclic groups to enable the decomposition of the source into 

1 We have restricted our attention to discrete memoryless sources. There has been a lot of research activity in the literature in 
the case of continuous-alphabet sources. Those works are not included in the reference list for conciseness. Please see the references 
in our earlier work on Gaussian sources [54] for a more complete list. 
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its constituent "digits" which are then encoded sequentially. Second, we show that, although group codes may 
not approach the Shannon rate-distortion function in a single source point-to-point setting, it is possible to 
construct non-trivial group codes which contain a code that approaches it. Using these two ideas, we provide an 
all-group-code solution to the problem and characterize an inner bound to the performance limit using single- 
letter information quantities. We also demonstrate the superiority of this approach over the conventional coding 
approach based on unstructured random codes for the case of reconstructing the modulo-2 sum of correlated 
binary sources with Hamming distortion. 

Special cases of the general problem of distributed source coding with joint distortion criterion have been 
studied before. The minimum rate at which a source X must be transmitted for the decoder to enable lossless 
reconstruction of a bivariate function F(X, Y) with perfect side information Y was determined in [7]. The 
case when the two communicators are allowed to exchange two messages was also considered. A two termi- 
nal interactive distributed source coding problem where the terminals exchange potentially infinite number of 
messages with the goal of reconstructing a function losslessly was studied in [8]. A similar problem of function 
computation from various sensor measurements in a wireless network was studied in [9] in the context of a 
packet collision model. 

Prior Work on Group Codes: Good codes over groups have been studied extensively in the literature when 
the order (size) of the group is a prime which enables the group to have a field structure. Such codes over 
Galois fields have been studied for the purpose of packing and covering (see [12,13] and the references therein). 
Two kinds of packing problems have received attention in the literature: a) combinatorial rigid packing and 
b) probabilistic soft packing, i.e., achieving the capacity of symmetric channels. Similarly, covering problems 
have been studied in two ways: a) combinatorial complete covering and (b) probabilistic almost covering, i.e., 
achieving the rate-distortion function of symmetric sources with Hamming distortion. Some of the salient 
features of these two approaches have been studied in [25]. In the following we give a sample of works in the 
direction of probabilistic packing and covering. Elias [1] showed that linear code achieve the capacity of binary 
symmetric channels. A reformulation of this result can be used to show [39] that linear codes can be used to 
losslessly compress any discrete source down to its entropy. Dobrushin [3] showed that linear codes achieve the 
random coding error exponent while Forney and Barg [18] showed that linear codes also achieve the expurgated 
error exponent. Further, these results have been shown to be true for almost all linear codes. Gallager [4] shows 
that binary linear codes succeeded by a nonlinear mapping can approach the capacity of any discrete memoryless 
channel. It follows from Goblick's work [2, 14, 15] on the covering radius of linear codes that linear codes can be 
used to achieve the rate distortion bound for binary sources with Hamming distortion. Blinovskii [16] derived 
upper and lower bounds on the covering radius of linear codes and also showed that almost all linear codes 
(satisfying rate constraints) are good source codes for binary sources with Hamming distortion. If the size of 
the finite field is sufficiently large, it was shown that in [17] that linear codes followed by a nonlinear mapping can 
achieve the rate distortion bound of a discrete memoryless source with arbitrary distortion measure. Wyner [40] 
derived an algebraic binning approach to provide a simple derivation of the Slcpian-Wolf [62] rate region for the 
case of correlated binary sources. Csiszar [53] showed the existence of universal linear encoders which attain 
the best known error exponents for the Slcpian-Wolf problem derived earlier using nonlinear codes. In [55,56], 
nested linear codes were used for approaching the Wyner-Ziv rate-distortion function for the case of doubly 
symmetric binary source and side information with Hamming distortion. Random structured codes have been 
used in other related multiterminal communication problems [19, 20, 67] to get performance that is superior 
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to that obtained by random unstructured codes. In [10], a coding scheme based on sparse matrices and ML 
decoding was presented that achieves the known rate regions for the Slepian-Wolf problem, Wyner-Ziv problem 
and the problem of lossless source coding with partial side information. 

Codes over general cyclic groups were first studied by Slepian [21] in the context of signal sets for the 
Gaussian channel. Forney [22] formalized the concept of geometrically uniform codes and showed that many 
known classes of good signal space codes were geometrically uniform. Biglieri and Elia [23] addressed the 
problem of existence of group codes for the Gaussian channel as defined by Slepian. Forney and Loeliger [24, 26] 
studied the state space representation of group codes and derived trellis representations which were used to build 
convolutional codes over abelian groups. An efficient algorithm for building such minimal trellises was presented 
in [27]. Loeliger [28] extended the concept of the Af-PSK signal set matched to the M-ary cyclic group to the 
case of matching general signal sets with arbitrary groups. Building codes over abelian groups with good error 
correcting properties was studied in [29]. The distance properties of group codes have also been extensively 
studied. In [30-32], bounds were derived on the minimum distance of group codes and it was also shown that 
codes built over nonabelian groups have asymptotically bad minimum distance behavior. Group codes have 
also been used to build LDPC codes with good distance properties [33]. The information theoretic performance 
limits of group codes when used as channel codes over symmetric channels was studied in [34] . Similar analysis 
for the case of turbo codes and geometrically uniform constellations was carried out in [35] . In [36] , Ahlswede 
established the achievable capacity using group codes for several classes of channels and showed that in general, 
group codes do not achieve the capacity of a general discrete memoryless channel. Sharper results were obtained 
for the group codes capacity and their upper bounds in [37,38]. 

The paper is organized as follows. In Section [2l we define the problem formally and present known results 
for the problem. In Section [3J we present an overview of the properties of groups in general and cyclic groups 
in particular that shall be used later on. We motivate our coding scheme in Section [U In Section we define 
the various concepts used in the rest of the paper. In Section [SI we present our coding scheme and present an 
achievable rate region for the problem defined in Section [2j Section [7] contains the various corollaries of the 
theorem presented in Section [5J These include achievable rates for lossless and lossy source coding while using 
group codes. We also present achievable rates using group codes for the problem of function reconstruction. 
Most of the proofs are presented in the appendix. In Section [5J we demonstrate the application of our coding 
theorem to various problems. We conclude the paper with some comments in Section [9[ 

A brief overview of the notation used in the paper is given below. Random variables are denoted by capital 
letters such as X, Y etc. The alphabet over which a discrete random variable X takes values will be indicated 
by X. The cardinality of a discrete set X is denoted by \X\. For a random variable X with distribution px('), 
the set of all n-length strongly e-typical sequences are denoted by A" (X) [57]. On most occasions, the subscript 
and superscript are omitted and their values should be clear from the context. For a pair of jointly distributed 
random variables X, Y with distribution px,y(', ■), the set of all n-length ^"-sequences jointly e-typical with a 
given x n sequence is denoted by the set A™(x n ). 

2 Problem Definition and Known Results 

Consider a pair of discrete random variables (X, Y) with joint distribution Pxy{',')- Let the alphabets of 
the random variables X and Y be X and y respectively. The source sequence (X n ,Y n ) is independent over 



4 



time and has the product distribution Pr((X n ,Y n ) = (x n ,y n )) = Y[7=i Vxy Hi)- We consider the following 
distributed source coding problem. The two components of the source are observed by two encoders which do 
not communicate with each other. Each encoder communicates a compressed version of its input through a 
noiseless channel to a joint decoder. The decoder is interested in reconstructing the sources with respect to a 
general fidelity criterion. Let Z denote the reconstruction alphabet, and the fidelity criterion is characterized by 
a mapping: d : X x y x Z — ► R+. We restrict our attention to additive distortion measures, i.e., the distortion 
among three n-length sequences x n , y n and z n is given by 

n 

d(x n ,y n ,z n ) = -^2 d ( Xi ^^- W 
n ■ i 

In this work, we will concentrate on the above distributed source coding problem (with one distortion 
constraint), and provide an information-theoretic inner bound to the optimal rate-distortion region. 

Definition 1. Given a discrete source with joint distribution Pxy{ x ,v) an d a distortion function d(-,-,-), a 
transmission system with parameters (n, 61,62, A) is defined by the set of mappings 

fi- X n —> {1, . . .,81}, f 2 :y n ^{l,...,6 2 } (2) 

g: {1,. .. A} x{l,. ..,&}->£" (3) 
such that the following constraint is satisfied. 

E(d(X n ,Y n ,g(MX n ),f 2 (Y n )))) < A. (4) 

Definition 2. We say that a tuple (R\, R 2 , D) is achievable i/Ve > 0, 3 for all sufficiently large n a transmission 
system with parameters (n, 61,82, A) such that 

-log6i<Ri + e fori = 1,2 A < D + e. (5) 

n 

The performance limit is given by the optimal rate- distortion region 1ZT> which is defined as the set of all 
achievable tuples (R\, R 2 , D). 

We remark that this problem formulation is very general. For example, defining the joint distortion measure 
d{X,Y,Z) as di(F(X,Y), Z) enables us to consider the problem of lossy reconstruction of a function of the 
sources as a special case. Though we only consider a single distortion measure in this paper, it is straightforward 
to extend the results that we present here for the case of multiple distortion criteria. This implies that the 
problem of reconstructing the sources subject to two independent distortion criteria (the Berger-Tung problem 
[46] ) can be subsumed in this formulation with multiple distortion criteria. The Slepian-Wolf [62] problem, the 
Wyner-Ziv problem [44], the Yeung-Berger problem [47] and the problem of coding with partial side information 
[42,63] can also be subsumed by this formulation since they all are special cases of the Berger-Tung problem. 
The problem of remote distributed source coding [6,43], where the encoders observe the sources through noisy 
channels, can also be subsumed in this formulation using the techniques of [48,49]. We shall see that our 
coding theorem has implications on the tightness of the Berger-Tung inner bound [46]. The two-user function 
computation problem of lossy reconstruction of Z = F(X, Y) can also be viewed as a special case of three-user 
Berger-Tung problem of encoding the correlated sources {X, Y, Z) with three independent distortion criteria, 
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where the rate of the third encoder is set to zero and the distortions of the first two sources are set to their 
maximum values. We shall see in Section 18.21 that for this problem, our rate region indeed yields points outside 
the Berger-Tung rate region thus demonstrating that the Berger-Tung inner bound is not tight for the case of 
three or more sources. 

An achievable rate region for the problem defined in Definitions [1] and [5] can be obtained based on the Berger- 
Tung coding scheme [46] as follows. Let V denote the family of pair of conditional probabilities (Pu\x> Pv\y) 
defined on X x U and y x V, where U and V are finite sets. For any (Pjj\x, Pv\y) £ P, let the induced joint 
distribution be Pxyuv = PxyP\j\xPv\y ■ U, V play the role of auxiliary random variables. Define G : U x V — » Z 
as that function of U, V that gives the optimal reconstruction Z with respect to the distortion measure d(-, ■, ■). 
With these definitions, an achievable rate region for this problem is presented below. 

Fact 1. For a given source {X, Y) and distortion d(-, •, •) define the region 1ZT>bt as 

TZV BT = (J Ir 1 >I(X;U\V),R 2 >I(Y;V\U),R 1 +R 2 >I(XY;UV), 

(Pu\x,Pv\v)eV I 

D >Ed(X,Y,G(U,V))^ (6) 
Then any (Ri, R 2l D) £ 7ZT> BT is achievable where * denotes convex closur^. 

Proof: Follows from the analysis of the Berger-Tung problem [46] in a straightforward way. □ 



3 Groups - An Introduction 

In this section, we present an overview of some properties of groups that are used later. We refer the reader 
to [59] for more details. It is assumed that the reader has some basic familiarity with the concept of groups. We 
shall deal exclusively with abelian groups and hence the additive notation will be used for the group operation. 
The group operation of the group G is denoted by +<?. Similarly, the identity element of group G is denoted 
by e<3. The additive inverse of a G G is denoted by —a. The subscripts are omitted when the group in question 
is clear from the context. A subset H of a group G is called a subgroup if H is a group by itself under the 
same group operation +g- This is denoted by H < G. The direct sum of two groups G\ and G 2 is denoted by 
Gi ®G 2 . The direct sum of a group G with itself n times is denoted by G n . 

An important tool in studying the structure of groups is the concept of group homomorphisms. 

Definition 3. Let G, H be groups. A function (f>: G — > H is called a homomorphism if for any a,b G G 

Cb(a+ G b)=0(a)+ H (b(b). (7) 
A bijective homomorphism is called an isomorphism. If G and H are isomorphic, it is denoted as G = H. 

A homomorphism (/>(•) has the following properties: 0(ec) = and a) = —0(a). The kernel ker(</>) of 
a homomorphism is defined as ker(</>) = {x 6 G: <fi(x) — e#}. An important property of homomorphisms is 

2 The cardinalities of U and V can be bounded using Caratheodary theorem [57]. 



(i 



that they preserve the subgroup structure. Let <\>: G —> H be a homomorphism. Let A < G and B < H. Then 
(f>^ 1 (B) < G and <fi(A) < B. In particular, taking B — {e#}, we get that ker(</>) < G. 

One can define a congruence result analogous to number theory using subgroups of a group. Let H < G. 
Consider the set Ha — {h + a : h G H}. The members of this set form an equivalence class called the right coset 
of H in G with a as the coset leader. The left coset of H in G is similarly defined. Since we deal exclusively 
with abelian groups, we shall not distinguish cosets as being left or right. All cosets are of the same size as H 
and two different cosets are either distinct or identical. Thus, the set of all distinct cosets of H in G form a 
partition of G. These properties shall be used in our coding scheme. 

It is known that a finite cyclic group of order n is isomorphic to the group Z„ which is the set of integers 
{0, . . . ,n — 1} with the group operation as addition modulo-n. A cyclic group whose order is the power of a 
prime is called a primary cyclic group. The following fact demonstrates the role of primary cyclic groups as the 
building blocks of all finite abelian groups. 

Fact 2. Let G be a finite abelian group of order n > 1 and let the unique factorization of n into distinct prime 
powers be n = Y\i=iPT ■ Then, 

G = Ai® A 2 - ■ ■ ® A k where \AA = p-' (8) 

Further, for each Ai,l < i < k with \Ai\ — pi' , we have 

Ai = Z p p © Z p H 2 • • • © Z p K t (9) 

where h\ > hi ■ ■ • > ht and X)*=i hj ~ e i- This decomposition of Ai into direct sum of primary cyclic groups is 
called the invariant factor decomposition of Ai . Putting equations and (0j together, we get a decomposition 
of an arbitrary abelian group G into a direct sum of possibly repeated primary cyclic groups. Further, this 
decomposition of G is unique,i.e., if G = B\ © B 2 ■ ■ ■ B m with \B.- t \ = pp for all i, then Bi = Ai and Bi and Ai 
have the same invariant factors. 

Proof: See [59], Section 5.2, Theorem 5. □ 

For example, Fact [2] implies that any abelian group of order 8 is isomorphic to either Zg or Z4 © Z2 or to 
Z2 © Z2 © Z2 where © denotes the direct sum of groups. Thus, we first consider the coding theorems only for 
the primary cyclic groups Z pr . Results obtained for such groups are then extended to hold for arbitrary abelian 
groups through this decomposition. Suppose G has a decomposition G = Z p «i © • • • © Z p =r where p\ > • • • > p r 
are primes. A random variable X taking values in G can be thought of as a vector valued random variable 
X = (Xi, . . . , X r ) with Xi taking values in the cyclic group Z p =i , 1 < i < r. Xi are called the digits of X. 

We now present some properties of primary cyclic groups that we shall use in our proofs. The group Z m is a 
commutative ring with the addition operation being addition modulo-m and the multiplication operation being 
multiplication modulo-m. This multiplicative structure is also exploited in the proofs. The group operation in 
Z 7 ^ is denoted by u™ + u 2 l - Addition of w™ with itself k times is denoted by ku\. The multiplication operation 
between elements x and y of the underlying ring Z m is denoted by xy. We shall say that x € Z m is invertible 
if there exists y G Z m such that xy — 1 where 1 is the multiplicative identity of Z m . The multiplicative inverse 
of x G Z m , if it exists, is denoted by a: -1 . The additive inverse of it™ G Z^ which always exists is denoted by 
— u™. The group operation in the group Z m is often explicitly denoted by ffi m . 
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We shall build our codebooks as kernels of homomorphisms from Z" r to Wl T . Justification for restricting the 
domain of our homomorphisms to Z™ r comes from the decomposition result of Fact [2j The reason for restricting 
the image of the homomorphisms to J]^ r shall be made clear later on (see the proof of Lemma 2]) . We need the 
following lemma on the structure of homomorphisms from Z™ r to Z^ r . 

Fact 3. Let Rom(fE^ r , Z^ r ) be the set of all homomorphisms from the group Z™ r to Z^ r and M(k,n,Z pr ) be the 
set of all k x n matrices whose elements take values from the group Z pr . Then, there exists a bisection between 
HomCZpr ,Zpr) and M(k,n,Z pr ) given by the invertible mapping f: i/om(Z" r , Z pr ) — > M(k,n, Z pr ) defined as 
f(4>) = such that <p(x n ) = $ • x" for all x n 6 Z™ r . Here, the multiplication and addition operations involved 
in the matrix multiplication are carried out modulo -p r . 

Proof: See [60], Section VI. □ 

4 Motivation of the Coding Scheme 

In this section, we present a sketch of the ideas involved in our coding scheme by demonstrating them for the 
simple case when the sources are binary. The emphasis in this section is on providing an overview of the main 
ideas and the exposition is kept informal. Formal definitions and theorems follow in subsequent sections. We 
first review the linear coding strategy of [39] to reconstruct losslessly the modulo-2 sum of Z = X ©2 Y of the 
binary sources X and Y. We then demonstrate that the Slepian-Wolf problem can be solved by a similar coding 
strategy. We generalize this coding strategy for the case when the fidelity criterion is such that the decoder 
needs to losslessly reconstruct a function F(X, Y) of the sources. This shall motivate the problem of building 
"good" channel codes over abelian groups. We then turn our attention to the lossy version of the problem 
where the sources X and Y are quantized to U and V respectively first. For this purpose, we need to build 
"good" source codes over abelian groups. Then, encoding is done in such a way that the decoder can reconstruct 
G(U, V) which is the optimal reconstruction of the sources with respect to the fidelity criterion •, •) given 
U, V. This shall necessitate the need for "good" nested group codes where the coarse code is a good channel 
code and the fine code is a good source code. These concepts shall be made precise later on in Sections and [5] 

4.1 Lossless Reconstruction of the Modulo-2 Sum of the Sources 

This problem was studied in [39] where an ingenious coding scheme involving linear codes was presented. This 
coding scheme can be understood as follows. It is well known [40] that linear codes can be used to losslessly 
compress a source down to its entropy. Formally, for any binary memory less source Z with distribution pz(z) 
and any e > 0, there exists a k x n binary matrix A with — < H(Z) + e and a function ip such that 

P{iP{Az n ) ^ z 11 ) < e (10) 

for all sufficiently large n. Let Z = X (B2 Y be the modulo-2 sum of the binary sources X and Y. Let the 
matrix A satisfy equation (jTUJ) . The encoders of X and Y transmit Si = Ax n and S2 — Ay n respectively at rates 
(H(Z),H(Z)). The decoder, upon receiving si ands2, computes ip(si (B2 S2) = tp(Ax n ®2Ay n ) =%p(Az n ). Since 
the A matrix was chosen in accordance with equation (JTUJ) , the decoder output equals z n with high probability. 
Thus, the rate pair (H(Z), H(Z)) is achievable. If the source statistics is such that H (Z) > H(X), then clearly 
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it is better to compress X at a rate H(X). Thus, the Korner-Marton coding scheme achieves the rate pair 
(Ri,R 2 ) with R x > mm{H(X),H{Z)} and R 2 > mm{H(Y), H(Z)}. This coding strategy shall be referred to 
as the Korner-Marton coding scheme from here on. 

The crucial part played by linear codes in this coding scheme is noteworthy. Had there been a centralized 
encoder with access to x n and y n , the coding scheme would be to compute z n — x n @2U n first and then compress 
it using any method known to achieve the entropy bound. Because the encoding is linear, it enables the decoder 
to use the distributive nature of the linear code over the modulo-2 operation to compute s\ ©2 S2 = Az n . Thus, 
from the decoder's perspective, there is no distinction between this distributed coding scheme and a centralized 
scheme involving a linear code. Also, in contrast to the usual norm in information theory, there is no other 
known coding scheme that approaches the performance of this linear coding scheme. 

More generally, in the case of a prime q, a sum rate of 2H{X ® q Y) can be achieved [58] for the reconstruction 
of the sum of the two g-ary sources Z = X@ q Y in any prime field 1 q . Abstractly, the Korner-Marton scheme can 
be thought of as a structured coding scheme with codes built over groups that enable the decoder to reconstruct 
the group operation losslessly. It turns out that extending the scheme would involve building "good" channel 
codes over arbitrary abelian groups. It is known (see Fact [5]) that primary cyclic groups Z p r are the building 
blocks of all abelian groups and hence it suffices to build "good" channel codes over the cyclic groups Z p .- . 

4.2 Lossless Reconstruction of the Sources 

The classical result of Slepian and Wolf [62] states that it is possible to reconstruct the sources X and Y 
noiselessly at the decoder with a sum rate of R\ + R2 — H(X,Y). As was shown in [53], the Slepian- Wolf 
bound is achievable using linear codes. Here, we present an interpretation of this linear coding scheme and 
connect it to the one in the previous subsection. We begin by making the observation that reconstructing the 
function Z = (X, Y) for binary sources can be thought of as reconstructing a linear function in the field Z2 ©Z2. 
This equivalence is demonstrated below. Let the elements of Z2 © Z2 be {00,01, 10, 11}. Denote the addition 
operation of Z 2 © Z 2 by ®k- 
Define the mappings 



Clearly, reconstructing (X, Y) losslessly is equivalent to reconstructing the function Z = X (Bk Y losslessly. 
The next observation is that elements in Z2©^2 can be represented as two dimensional vectors whose components 
are in Z2. Further, addition in Z2 © Z2 is simply vector addition with the components of the vector added in 
Z2. Let the first and second bits of X be denoted by X\ and X 2 respectively. The same notation holds for Y 
and Z as well. Then, we have the decomposition of the vector function Z as Zi — Xi ©2 Yi for i = 1,2. 

Encoding the vector function Z directly using the Korner-Marton coding scheme would entail a sum rate of 
Ri+ R 2 = min{H(X,Y),H(X)} + min{H(X, Y), H(Y)} = H{X) + H(Y) which is more than the sum rate 
dictated by the Slepian- Wolf bound. Instead, we encode the scalar components of the function Z sequentially 
using the Korner-Marton scheme. Suppose the first digit plane Z\ is encoded first. Assuming that it gets 




(11) 



(12) 
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decoded correctly at the decoder, it is available as side information for the encoding of the second digit plane 
Z 2 . Clearly, the Korner-Marton scheme can be used to encode the first digit plane Z x , The rate pair (i?n, R 2X ) 
achieved by the scheme is given by 

Ru > mm{H{Zi), H(X X )} = H{X X ) = (13) 
R 21 > min{H(Z x ), H(Y X )} = H{Z X ) (14) 

It is straightforward to extend the Korner-Marton coding scheme to the case where decoder has available to 
it some side information. Since Z\ is available as side information at the decoder, the rates needed to encode 
the second digit plane Z 2 are 

R 12 > mm{H(Z 2 \ Z X ),H(X 2 \ Z x )} = H(Z 2 \ Z x ) (15) 
R 22 > min{H(Z 2 \ Z X ),H(Y 2 \ Z x )} = H(Y 2 \ Z x ) = (16) 

Thus, the overall rate pair needed to reconstruct the sources losslessly is 

R x = R xx + R X2 > H(Z 2 I Z x ) = H(X 2 \ Y x ) (17) 
R 2 =R 2X + R 22 >H{Z X )=H(Y X ), (18) 

The sum rate for this scheme is R\ + R 2 = H(X 2 ,Y X ) = H (X, Y) thus equaling the Slepian-Wolf bound. 

4.3 Lossless Reconstruction of an Arbitrary Function F(X, Y) 

While there are more straightforward ways of achieving the Slepian-Wolf bound than the method outlined 
in Section 14.21 our encoding scheme has the advantage of putting the Korner-Marton coding scheme and the 
Slepian-Wolf coding scheme under the same framework. The ideas used in these two examples can be abstracted 
and generalized for the problem when the decoder needs to losslessly reconstruct some function F(X, Y) in order 
to satisfy the fidelity criterion. 

Let us assume that the cardinality of X and Y are respectively a and (3. The steps involved in such an 
encoding scheme can be described as follows. We first represent the function as equivalent to the group operation 
in some abelian group A. This is referred to as "embedding" the function in A. This abelian group is then 
decomposed into its constituent cyclic groups and the embedded function is sequentially encoded using the 
Korner-Marton scheme outlined in Section I4TT1 Encoding is done keeping in mind that, to decode a digit, the 
decoder has as available side information all previously decoded digits. 

It suffices to restrict attention to abelian groups A such that \Z\ < \A\ < a(3. Clearly, if the function 
Fi(X, Y) = (X, Y) can be embedded in a certain abelian group, then any function F(X, Y) can be reconstructed 
in that abelian group. This is because the decoder can proceed by reconstructing the sources {X, Y) and then 
computing the function F(X, Y). It can be shown (see Appendix [E| that the function F X {X, Y) = (X, Y) can 
be reconstructed in the group Z a © which is of size a/3. Clearly, \A\ > \Z\ is a necessary condition for the 
reconstruction of Z = F(X,Y). 

4.4 Lossy Reconstruction 

We now turn our attention to the case when the decoder wishes to obtain a reconstruction Z with respect to 
a fidelity criterion. The coding strategy is as follows: Quantize the sources X and Y to auxiliary variables U 
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and V. Given the quantized sources U and V, let G(U, V) be the optimal reconstruction with respect to the 
distortion measure d(-, •,•). Reconstruct the function G(U,V) losslessly using the coding scheme outlined in 
Section [O] 

We shall use nested group codes to effect this quantization. Nested group codes arise naturally in the area of 
distributed source coding and require that the fine code be a "good" source code and the coarse code be a "good" 
channel code for appropriate notions of goodness. We have already seen that to effect lossless compression, the 
channel code operates at the digit level, ft follows then that we must use a series of nested group codes, one 
for each digit, over appropriate cyclic groups. For instance, if the first digit of G(U, V) is over the cyclic group 
Z p «i , then we need nested group codes over Z p <=i that encode the sources X and Y to U\ and V\ respectively. 
The quantization operation is also carried out sequentially, i.e., the digits JJ% and V2 are encoded given the 
knowledge that either Z\ or (Ui,Vi) is available at the decoder and so on. The existence of "good" nested 
group codes over arbitrary cyclic groups is shown later. 

The steps involved in the overall coding scheme can be detailed as follows: 

• Let U, V be discrete random variables over the alphabet U, V respectively. Further suppose that \U\ = 
a, I V| = j3. Choose the joint density Px,y,u,v — Px,yPu\xPv\y satisfying the Markov chain U — X — Y — V . 

• Let G(U, V) be the optimal reconstruction function with respect to d(-, ■, ■) given U, V. 

• Embed the function G(U, V) in an abelian group A, \Q\ < \A\ < a(3. 

• Decompose G(U, V) into its constituent digit planes. Fix the order in which the digit planes are to be 
sequentially encoded. 

• Suppose the digit plane is the cyclic group 1> p " b . Quantize the sources (X n ,Y n ) into digits (Ub,Vb) 
using the digits already available at the decoder as side information. The details of the quantization 
procedure are detailed later. 

• Encode Zf, = Uj, ® p e b Vb using group codes. 

5 Definitions 

When a random variable X takes value over the group Z pr , we need to ensure that it doesn't just take values 
in some proper subgroup of Z r . This leads us to the concept of a non-redundant distribution over a group. 

Definition 4. A random variable X with X = % r or its distribution Px is said to be non-redundant if 
Px{x) > for at least one symbol x € Z pr \pZ pr . 

It follows from this definition that x n E A"(X) contains at least one x € 7L pT \p1 pT if X is non-redundant. 
Such sequences are called non-redundant sequences. A redundant random variable taking values over 1 pT can be 
made non-redundant by a suitable relabeling of the symbols. Also, note that a redundant random variable over 
7j pr is non-redundant when viewed as taking values over 1 p T-i for some < i < r. Our coding scheme involves 
good nested group codes for source and channel coding and the notion of embedding the optimal reconstruction 
function in a suitable abelian group. These concepts are made precise in the following series of definitions. 
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Definition 5. A bivariate function G: U x V — > Q is said to be embeddable in an abelian group A with respect 
to the distribution puv{u,v) on IA x V if there exists injective functions : U — > A, Sy 4 ' : V — > A and a 
surjective function : yl — > Q such that 

S { c ) (S i u A} (u)+ A S i v A) (v))=G(u,v) V(u,«) E W x V with Puv (u,v) > (19) 

// G(J7, V) is indeed embeddable in the abelian group A, it is denoted as G{U, V) C A with respect to the 
distribution puv{u,v). Define the mapped random variables U = S{f~\u) and V = Sy (V). Their dependence 
on A is suppressed and the group in question will be clear from the context. 

Suppose the function G(U,V) C A with respect to puv- We encode the function G(U,V) sequentially by 
treating the sources as vector valued over the cyclic groups whose direct sum is isomorphic to A. This alternative 
representation of the sources is made precise in the following definition. 

Definition 6. Suppose the function G(U,V) C A with respect to puv ■ Let A be isomorphic to ®i_-yZ p H where 
Pi < ■ • ■ < Pk are primes and ej are positive integers. Then, it follows from Fact [H that there exists a bisection 
Sa ■ A — ► Z p «i x . . . Z p i fc . Let U = Sa(U), V — Sa(V). Let U — (Ui, . . . , Uk) be the vector representation of U. 
The random variables Ui are called the digits of U. A similar decomposition holds for V . Define Z — [Z\ , . . . , Z/.) 
where Z { = U l V l . It follows that S A \Z) = U + A V. 

Our encoding operation proceeds thus: we reconstruct the function G(U, V) by first embedding it in some 
abelian group A and then reconstructing U +a V which we accomplish sequentially by reconstructing Ui © p <=; Vi 
one digit at a time. While reconstructing the ith digit, the decoder has as side information the previously 
reconstructed (i — 1) digits. This digit decomposition approach requires that we build codes over the primary 
cyclic groups 1i pr which are "good" for various coding purposes. We define the concepts of group codes and 
what it means for group codes to be "good" in the following series of definitions. 

Definition 7. Let A be a finite abelian group. A group code C of blocklength n over the group A is a subset 
of A n which is closed under the group addition operation, i.e., C C A n is such that i/c",C2 E C, then so does 

Cl +a* C§. 

Recall that the kernel ker(0) of a homomorphism (f>: A n — * A k is a subgroup of A n . We use this fact to build 
group codes. As mentioned earlier, we build codes over the primary cyclic group 1 pT . In this case, every group 
code C C ZpV has associated with it a k x n matrix H with entries in Z r which completely defines the group 
code as 

C = {x n E ZpV: Hx n = fc }. (20) 

Here, the multiplication and addition are carried out modulo-p r . H is called the parity-check matrix of the 
code C. We employ nested group codes in our coding scheme. In distributed source coding problems, we often 
need one of the components of a nested code to be a good source code while the other one to be a good channel 
code. We shall now define nested group codes and the notions of "goodness" used to classify a group code as a 
good source or channel code. 

Definition 8. A nested group code {C\,C2) is a pair of group codes such that every codeword in the codebook 
C2 is also a codeword in C±, i.e., C2 <C\. Their associated parity check matrices are the k\ X n matrix Hi and 
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the ki x n matrix Hi. They are related to each other as H\ — J ■ Hi for some k\ x ki matrix J. One way to 
enforce this relation between Hi and Hi would be to let 



H 2 



H x 
AH 



(21) 



where AH is a (k 2 — k\) x n matrix over Z r . 



The code C\ is called the fine group code while Ci is called the coarse group code. When nested group codes 
are used in distributed source coding, typically the coset leaders of Ci in C\ are employed as codewords. In such 
a case, the rate of the nested group code would be n~ 1 (ki — k\) logp r bits. 

We define the notion of "goodness" associated with a group code below. To be precise, these notions are 
defined for a family of group codes indexed by the blocklcngth n. However, for the sake of notational convenience, 
this indexing is not made explicit. 

Definition 9. Let Pxu be a distribution over X x U such that the marginal Pjj is a non-redundant distribution 
over Z pr for some prime power p r . For a given group code C over hi and a given e > 0, let the set A e {C) be 
defined as 

A € (C) = {x n : 3u n e C such that (x n , u n ) e A™(X, U)}. (22) 
The group code C overhi is called a good source code for the triple (X,U,Pxu) if we have Ve > 0, 

P%(A e (C))>l-e (23) 

for all sufficiently large n. 

Note that, a group code which is a good source code in this sense may not be a good source code in the 
usual Shannon sense. Rather, such a group code contains a subset which is a good source code in the Shannon 
sense for the source Px with forward test channel Pjj\x- 

Definition 10. Let Pzs be a distribution over Z xS such that the marginal Pz is a non-redundant distribution 
over Is p r for some prime power p r . For a given group code C over Z and a given e > 0, define the set B e (C) as 
follows: 

B e {C) = {(z n , s n ) : 3~z n such that (F\ s n ) E A^{Z, S) and H~z n = Hz n }. (24) 

Here, H is the k(n) x n parity check matrix associated with the group code C. The group code C is called a good 
channel code for the triple (Z,S, Pzs) if we have Ve > 7 

FSs(B e {C)) < e (25) 

for all sufficiently large n. Associated with such a good group channel code would be a decoding function ip : 
7L k r x S n -> U v \ such that 

P(iP(Hz n ,s n ) = z n )>l-e. (26) 

Note that, as before, a group code which is a good channel code in this sense may not a good channel code 
in the usual Shannon sense. Rather, every coset of such a group code contains a subset which is a good channel 
code in the Shannon sense for the channel Ps\z with input distribution Pz- This interpretation is valid only 
when S is a non-trivial random variable. 
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Lemma 1. For any triple (Z,S, Pzs) °f two finite sets and a distribution, with \Z\ = p r a prime power and 
Pz non- redundant, there exists a sequence of group codes C that is a good channel code for the triple {Z, S, Pzs) 
such that the dimensions of their associated k(n) x n parity check matrices satisfy 

lira M^>logp r = max ( —) {H{Z\S) - H{[Z] t \S)) (27) 



n— »oo n 0<i<r \r 

where \Z\ is a random variable taking values over the set of all distinct cosets of p ir L pr in Z^r. For example, if 
Z = 7js, then \Z\i is a A-ary random variable with symbol probabilities (pz(0) +Pz(4)), (pz(l) +pz(5)), (pz(2) + 
pz{6)) and (p z {3)+p z {7)). 

Proof: See Appendix \X[ □ 

Note that [Z]o is a constant and [Z] T = Z. When building codes over groups, each proper subgroup of the 
group contributes a term to the maximization in equation (|27[) . Since the smaller the right hand side of equation 
(|27jl . the better the channel code is, we incur a penalty by building codes over groups with large number of 
subgroups. 

Lemma 2. For any triple {X,U,Pxu) of two finite sets and a distribution, with \U\ — p r a prime power and 
Pu non-redundant, there exists a sequence of group codes C that is a good source code for the triple (X ,U,Pxu) 
such that the dimensions of their associated k(n) x n parity check matrices satisfy 

Km Ml^iogp 1 - = min( H(U\X),r\H(U\X) -logp r - 1 | + ) (28) 
n— *oo n 

where \x\ + = max(i, 0). 

Proof: See Appendix [51 □ 

Putting r = 1 in equations ([2"T[) and (|28p, we get the performance obtainable while using linear codes built 
over Galois fields. 

Lemma 3. Let X, Y, S, U, V be five random variables where U and V take value over the group Z pr for some 
prime power p r . Let Z = U ® p r V . Let U — > X — > Y — > V form a Markov chain, and let S — > (X, Y) — > (U, V) 
form a Markov chain. From the Markov chains, it follows that H(U\X) < H(Z\S), H(V\Y) < H(Z\S). Without 
loss of generality, let H(U\X) < H(V\Y) < H{Z\S). Then, there exists a pair of nested group codes (Cn,C2) 
and (Ci2,C2) such that 

• Cn is a good group source code for the triple {X,U, Pxu) with 

lim kl1 ^ log;/ = mm(H(U\X),r\H(U\X) - logp^ 1 ^) (29) 

n — >oo n 



C12 is a good group source code for the triple (y, V, Pyv) with 



li m M^iogp r =mm{H(V\Y),r\H(V\Y) -logp^ 1 ^) (30) 

n— >oo n 



C2 is a good group channel code for the triple (Z,S,Pzs) with 



li m M^l l ogp r = max ( — ) {H{Z\S) - H{[Z] t \S)) (31) 



n— >CX3 n 0<i<r \ r — I 
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Proof: See Appendix [Cl 



□ 



Note that while choosing the codebooks C\\, Cyi and C2, the perturbation parameters e in Definitions [51 andfTUl 
need to be chosen appropriately relative to each other so that the n-length sequences (X n , Y n , S n , U n , V n , Z n ) 
are jointly typical with high probability. Due to the Markov chains U — » X — ► Y — » V and S — > (X, F) — * ([/, V), 
it follows from Markov lemma [51] that if (X n , Y n , S n ) is generated according to Pxys and if U n is generated 
jointly typical with X n and V" is generated jointly typical with Y n , then (X n , Y n , S n , U n , V n , Z n ) is jointly 
strongly typical (for an appropriate choice of e) with high probability. 



6 The Coding Theorem 

We are given discrete random variables X and Y which are jointly distributed according to Pxy ■ Let V denote 
the family of pair of conditional probabilities {Pxj\Xi Pv\y) defined on X x U and y x V, where U and V are finite 
sets, \U\ = a, |V| = (3. For any (Pu\x>Pv\y) 6 P ■> ^ the induced joint distribution be Pxyuv — PxyPu\xPv\y- 
U, V play the role of auxiliary random variables. Define G: U x V — * Z as that function of U, V that gives the 
optimal reconstruction Z with respect to the distortion measure d(-, ■, •). Let Q denote the image of G(U,V). 
Let T = {A: A is abelian, \Q\ < \A\ < af3 : G(U, V) <Z A with respect to Pjjv}- It is shown in Appendix [El that 
the set T is non-empty, i.e., there always exists an abelian group A G T in which any function G(U,V) can 
be embedded. For any A G T, let A be isomorphic to ®!? =1 Z p H. Let U = S A (Sif ] (U)) and V = S A (S^\v)) 
where the mappings are as defined in Definitions O and [51 Define Z = (Z\, . . . , Zk) where Zi = Ui® % and the 
addition is done in the group to which the digits Ui , Vi belong. Assume without loss of generality that the digits 
Ui, Vi, Zi, 1 < i < k are all non-redundant. If they are not, they can be made so by suitable relabeling of the 
symbols. Recall the definition of [Z]i from Lemma[TJ The encoding operation of the X and Y encoders proceed 
in k steps with each step producing one digit of U and V respectively. Let tta- {1, . . . , fc} — > {1, . . . , k} be a 
permutation. The permutation tta can be thought of as determining the order in which the digits get encoded 
and decoded. Let the set 11,4(6), 1 < b < k be defined as 11^(6) = {I: tta(I) < b). The set 11^(6) contains the 
indices of all the digits that get encoded before the 6th stage. At the 6th stage, let the digits U 7TA m , V 7TA m take 
values over the group Z££. With these definitions, an achievable rate region for the problem is presented below. 

Theorem 1. For a given source (X, Y), define the region lZT>i n as 

( k k 

TZV m 4 (J J (R u R 2 , D) : R ± > ^ min (r$ , R$ ) , R 2 > £ min (b£} , R%j} ) (32) 



r U\X <^V\Y>' 
A£T,TT A 



6=1 b=l 



D > Ed(X, Y, G(U, V)) } (33) 



whe 



^16 > 



(b)}i\Zu A (b)) 

min(i/(C/^ (&) |X,C/n A (6)),r6(| J ff(t>^ (fc) I A,?7n A (6))-logpr _1 | + ))] (34) 
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and 



R i2) > 



max 



0<i<r b \r b — I 



n 



(•ff(^7r A (6) I z n A (b)) - H([U 7rA ( b )]i I Z UA ^ b ))j 



min (H{U„ m | X,£ nA(6) ), r 6 (|ff(^ (6) | X, U^ m ) - log^" 1 ^))] (35) 

The quantities and R^ b are similarly defined with (X, U) replaced by (Y, V). Then any (i?i, R 2 , D) S TLT)* n 
is achievable where * denotes convex closure. 

Proof: Since the encoders don't communicate with each other, we impose the Markov chain V — Y — X — U 
on the joint distribution Pxyuv- The family V contains all distributions that satisfy this Markov chain. Fix 
such a joint distribution. Fix A £ T and the permutation it a' {1, ■ ■ • , k} —* {1, . . . , k}. The encoding proceeds 
in k stages with the &th stage encoding the digits U 7rA ^ b ) 1 V„ A n,) m order to produce the digit Z 7TA ^y For this, 
the decoder has side information Zn A ( b ) . 

Let U 7TA , V WA take values over the group 17* b . The encoders have two encoding options available at the 
&th stage. They can either encode the digits U 7TA (b) and V-K A (b) directly or encode in such a way that the decoder 
is able to reconstruct Z VA i b \ directly. We present a coding scheme to achieve the latter first. 

We shall use a pair of nested group codes (Cub,C 2 b) and (Ci 2 b,C 2 b) to encode Z^ A t b y Let the corresponding 
parity check matrices of these codes be Hub, H\ib and H 2b respectively. Let the dimensionality of these matrices 
be k\\\, x n, kyib x n and k2b x n respectively. These codebooks are all over the group Z£*. We need Cub to be a 
good source code for the triple {X x Un A (b), £Y WA ( & ), Pxu u {b) u (b) )> ^ 12f > to De a g°°d source code for the triple 
CVxVn^fe), V WA ( 6 ), P YV UAm % Am ) andC 2fc to be a good channel code for the triple (Z nA ( b ), Zn A (b), p z^ A(b) z UA(b) )- 

The encoding scheme used by the X-encoder to encode the &th digit, 1 < b < k is detailed below. The 
X-encoder looks for a typical sequence G Cub such that it is jointly typical with the source sequence X™ 

and the previous encoder output digits U^ A ( b y If it finds at least one such sequence, it chooses one of these 
sequences and transmits the syndrome Sxb = #26^™ a (m to the decoder. If it finds no such sequence, it declares 
an encoding error. The operation of the K-encoder is similar. 

Let tpb(',') be the decoder corresponding to the good channel code C 2 b- The decoder action is described by 
the following series of equations. The decoder receives the syndromes Sxb and Syb- 



Z. 



7ta(&) = (Sx b ® Pb r b Sy b , Zn A{b )) 



V>& (H 2b U. 

i> b (H 2b ( 



7T A (b) ®Pb r >> #26K A (6)> z n A (b) 



U n A (b) ®Pb r » V K A (b) j > Z TL A (b) 



n ^T, A (b)^n A (b) 

= Z™ A (b) w ith high probability (36) 

where (a) follows from the fact that C 2 b is a good channel code for the triple (Z„ A ( b y Zn A (&), P z {b) z n (b) )- 

The rate expended by the X-encoder at the &th stage can be calculated as follows. Since Cub is a good 
source code for the triple (X x l^n A (b)T^Tr A (b)j P xu n (&)& we nave f rorn equation (|28p that the dimensions 

of the parity check matrix Hub satisfy 

^logp r b b < min (H{U VA{b) \ X,U„ Aib) ),r b (\H(U nA{b) | X, U UA{b) ) - log^" 1 ^)) - e a (37) 
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Since C26 is a good channel code for the triple (Z 7TA ^ , Zn A (b) i P% ^ z n ( , ) ) > the dimensions of the parity check 
matrix H2b satisfy 

^log# > max (-^:) (#(4 A (b) I ^ru(b)) - ^([4 A (6)]«|^( 6 ))) + £2 (38) 
The rate of the nested group code in bits would be i?i = n~ l {k2b — &iib) logp^ 6 . Therefore, 



-""lb — 



max ' 

0<i<r 6 V 7*& — z 



H ( Z ir A (b) I ^n A (b)) _ -^([• Z ^7r J 4(6)]i|^n A (b)) 



min (if (# WA(6) |X,^n A(6) ),r 6 (| J ff(C/ 7rA(fe) | X, ^(b)) - log^'V))] + e 1 + e 2 (39) 

The other option that the encoders have is to directly encode the digits U 7VA ^ and V 7VA ^y This can also be ac- 
complished using nested group codes as follows. The X encoder uses the nested group code (Cub, C21&) such that 
the fine group code Ciu is a good source code for the triple {X x Un A (b)i Mn A (b) > ^xu u {b) u {b) ) an< ^ ^ 21b * s a g°°d 
channel code for the triple (i4- A (b), -Zn A (&)> ^t/ (6) z n (6) )- The Y encoder uses the nested group code (Ci2b, C22b) 
such that the fine group code C i2 b is a good source code for the triple (y x Vji A (b)> V WA (a), -fyvh (6) i> (b) ) an d 
C22b is a good channel code for the triple ( V nA (p) , Zn A ^ , Py (b) z n The encoding operation is similar to 
that described earlier and it is easy to verify its correctness. 

The rate of this nested group code in bits would be R\ = n~ 1 (k2b — fciib) logp^ 6 . Therefore, 



0(2) 
K ib > 



max 



0<i<r b \Tb — I 



n, 



H{Un A (£) I z n A (b)) ~ -ff([#7r A (l>)]t I ^n A (b))^ 

- min (h (^ waW |X,r7 nA(6) ),rb(|ff(C/ 7rA( 6) | U nA (b)) ~ log^" 1 ^))] + e l + e 2 (40) 

Combining equations (|59")l and (|4H|) . we have proved Theorem [1] □ 

Remark 1: The design of the channel code used in the above derivation assumes that the side information 
available to the decoder at the 6th stage is Zjj A (b)- However, it is possible that at some stage 1 < i < k, the 
encoding was done in such a way that the decoder could decode (U 7rA (fj,V nA u)) and not just Z nA ^y Taking 
such considerations into account while designing the channel code for the 6th stage would lead to a possible 
improvement of the rate region in Theorem [1] 

Remark 2: In the above derivation, if the encoders choose to encode the sources U nA n,), V nA (b) directly instead 
of encoding the function Z 7rA (b), further rate gains are possible when one encoder encodes its source conditional 
on the other source in addition to the side information already available at the decoder. Such improvements 
are omitted for the sake of clarity of the expressions constituting the definition of the achievable rate region. 
Remark 3: The above coding theorem can be extended to the case of multiple distortion constraints in a 
straightforward fashion. 



7 Special cases 

In this section, we consider the various special cases of the rate region presented in Theorem [1] 
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7.1 Lossless Source Coding using Group Codes 

We start by demonstrating the achievable rates using codes over groups for the problem of lossless source coding. 
A good group channel code C for the triple (X, 0, Px) as defined in Definition [10] can be used to achieve lossless 
source coding of the source X. The source encoder outputs Hx n where H is the k x n parity check matrix of 
C. The decoder uses the associated decoding function ip(-, ■) to recover ip(Hx n ,0) = x n with high probability. 
From equation (|27[) . it follows that the dimensions of the parity check matrix satisfy 

- log^ > max (-!—) (H(X) - (41) 

Recognizing the term in the left as the rate of the coding scheme, we get the following corollary to Theorem [1] 

Corollary 1. Suppose X is a non redundant random variable over the group Z pr and the decoder wants to 
reconstruct X losslessly. Then, there exists a group based coding scheme that achieves the rate 

R > max ( — ^— ) (H(X) - H([X]i)) (42) 

0<i<r \r — i J 

Putting r = 1 in equation (]42[) reduces it to the well known result that linear codes over prime fields can 
compress a source down to its entropy. Thus, this achievable rate region using group codes can be strictly 
greater than Shannon entropy. A sufficient condition for the existence of group codes that attain the entropy 
bound is that 

H([X}i) > -H{X) for < i < r (43) 
r 

7.2 Lossy Source Coding using Group Codes 

We next consider the case of lossy point to point source coding using codes built over the group Z pr . . Consider 
a memoryless source X with distribution Px- The decoder attempts to reconstruct U that is within distortion 
D of X as specified by some additive distortion measure d: X xU — » R + . Suppose U takes its values from the 
group Z p ,.. A good group source code C for the triple (X,U, Pxu) as defined in Definition can be used to 
achieve lossy coding of the source X provided the joint distribution Pxu is such that E(d(X, U)) < D and U 
is non-redundant. The source encoder outputs u n £ C that is jointly typical with the source sequence x n . An 
encoding error is declared if no such u n is found. The decoder uses u n as its reconstruction of the source x n . 
From equation (|2"5|) , it follows that the dimensions of the parity check matrix associated with C satisfy 

-logp r < min( H(U\X),r\H(U\X) -logp r - 1 | + ). (44) 
n 

The rate of this encoding scheme is R = (l — -) \ogp r . Thus, we get the following corollary to Theorem[T] 

Corollary 2. Let X be a discrete memoryless source and U be the reconstruction alphabet. Suppose U = Z pr 
and the decoder wants to reconstruct the source to within distortion D as measured by the fidelity criterion d(-, •). 
Without loss of generality, assume that U is non-redundant. Then, there exists a group based coding scheme 
that achieves the rate 

R> min log/ - (min (H(U\X),r\H(U\X) - logp'" 1 ^)) . (45) 

p u\x 
Ed(X,U)<D 
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If U takes values in a general abelian group of order n that is not necessarily a primary cyclic group, then 
a decomposition based approach similar to the one used in the proof of Theorem [1] can be used. Suppose 
n = J\i = i Pi* is the prime factorization of n. Then, the group in which U takes values is isomorphic to ©|L 1 Z p fi 
where pi are not necessarily distinct primes. The random variable U can be decomposed into its constituent digits 
(U\, . . . , Uk) which can then be encoded sequentially. The achievable rate can be obtained in a straightforward 
way. A simplification occurs when U is treated as a subset of the group ©*L 1 Z Pi where pi are prime. In this 
case, the following rate-distortion bound is obtained. 

Corollary 3. Let X be a discrete memoryless source and U be the reconstruction alphabet. Letpx,...,pk be 
primes such that ]X=i£ , « — 1^1- Suppose the decoder wants to reconstruct the source to within distortion D as 
measured by the fidelity criterion d(-, •). Then, there exists a group based coding scheme that achieves the rate 

R> min ( V log Pi | - H(U\X) (46) 

Ed(X,U)<D ' 

Corollary [3] can be viewed as providing an achievable rate-distortion pair for lossy source coding using linear 
codes built over Galois fields. Note that it is possible to construct codebooks with rate R = H(U) — H(U\X) 
by choosing codewords independently and uniformly from the set A™(U). By imposing the group structure on 
the codebook, we incur a rate loss of (X}<=i ^°SPi ~~ H(U)) bits per sample. This rate loss is strictly positive 
unless the random variable U is uniformly distributed over ®\ = {l> Vi ■ 

7.3 Nested Linear Codes 

We specialize the rate region of Theorem [T] to the case when the nested group codes are built over cyclic groups 
of prime order, i.e., over Galois fields of prime order. In this case, group codes over Z pr reduce to the well known 
linear codes over prime fields. It was already shown in Sections 17.11 and 17.21 that Lemmas Q] and [2] imply that 
linear codes achieve the entropy bound and incur a rate loss while used in lossy source coding. In this section, 
we demonstrate the implications of Theorem [T] when specialized to the case of nested linear codes, i.e., when r 
is set to 1. 

7.3.1 Shannon Rate-Distortion Function 

We remark that Theorem [1] shows the existence of nested linear codes that can be used to approach the rate- 
distortion bound in the single-user setting for arbitrary discrete sources and arbitrary distortion measures. 

Corollary 4. Let X be a discrete memoryless source with distribution Px and let X be the reconstruction 
alphabet. Let the fidelity criterion be given by d: X x X — > R + . Then, there exists a nested linear code {Cx^C^) 
that achieves the rate- distortion bound 

R(D) = min I(X;X) (47) 

p - 

x\x 
Ed(X,X)<D 

Proof: Let the optimal forward test channel that achieves the bound be given by P x i x - Suppose q is a prime 
such that X C Z 9 and X is non- redundant. The rate bound, given by I(X;X) can be approached using a 
nested linear code (C\,C-2) built over the group Z q . Here C\ is a good source code for the triple (X, X,P XX ) 
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and C2 is a good channel code for the triple (X,S, Px S ) where S = {0} and S is a degenerate random variable 
with Ps(0) = 1. It follows from Lemmas [2] and [T] that the dimensions of the parity check matrices associated 
with C\ and C2 satisfy 

li m ^)-logq = H(X\X) (48) 

n— >oo Tl 

lim ^±\ ogq = H {X) (49) 

n — >oo fl 

Thus, the rate achieved by this scheme is given by n~ 1 (k2(n) — fci(n)) logq = I(X; X). □ 

This can be intuitively interpreted as follows. For a code to approach the optimal rate-distortion function, 
the "Voronoi" region (under an appropriate encoding rule) of most of the codewords should have a certain shape 
(say, shape A), and a high-probability set of codewords should be bounded in a region that has a certain shape 
(say, shape B). We choose C\ such that the "Voronoi" region (under the joint typicality encoding operation 
with respect to p x x ) of each codeword has shape A. C 2 is chosen such that its "Voronoi" region has shape 
B. Hence the set of "coset leaders" of C\ in C 2 forms a code that can approach the optimal rate-distortion 
function. This reminds us of a similar phenomenon first observed in the case of Gaussian sources with mean 
squared error criterion in [66] , where the performance of a quantizer is measured by so-called granular gain and 
boundary gain. Granular gain measures how closely the Voronoi regions of the codewords approach a sphere, 
and boundary gain measures how closely the boundary region approaches a sphere. 

7.3.2 Berger-Tung Rate Region 

We now show that Theorem [1] implies that nested linear codes built over prime fields can achieve the rate region 
of the Berger-Tung based coding scheme presented in Lemma [TJ 

Corollary 5. Suppose we have a pair of correlated discrete sources [X, Y) and the decoder is interested in 
reconstructing Z to within distortion D as measured by a fidelity criterion d: X x y x Z — ► R + . For this 
problem, an achievable rate region using nested linear codes is given by 

UV BT = U {{Ri,R 2 )-- Ri > I(X;U\Y), 

(Pu\x-P V \Y)ev 

R 2 > I{Y; V\X),Rx + R 2 > I(X; U) + I(Y; V) - I(U; V)} (50) 

where V is the family of all joint distributions Pxyuv that satisfy the Markov chain U — X — Y — V such that 
the distortion criterion Ed(X,Y, Z(U,V)) < D is met. Here Z(U,V) is the optimal reconstruction of Z with 
respect to the distortion criterion given U and V . 

Proof: We proceed by first reconstructing the function G(U, V) = (U, V) at the decoder and then computing 
the function Z(U, V). For ease of exposition, assume that U = V = 7L q for some prime q. If they are not, a 
decomposition based approach can be used and the proof is similar to the one presented below. Clearly, G(U, V) 
can be embedded in the abelian group A = Z 9 © Z 9 . The associated mappings are given by U = (U,0) and 
V = (0, V) where is the identity element in 7L q . Thus, Z x = U + = U and Z 2 = + V = V. Encoding is 
done in two stages. Let the permutation tta(') be the identity permutation. Substituting this into equations 
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(|59"]l and (gDD gives us 



R 11 >min{H(Z 1 ),H{U 1 )}-H(U 1 \X)=I(X;U 1 ) = I(X-U) 1 

R 21 >min{H(Z 1 ),H(V 1 )}-H(V 1 \Y)=0, 

R 12 > min{H(Z 2 \ Z 1 ), H(U 2 \ Ux)} - H(U 2 \ X, U x ) = 0, 

R 22 > mm{H(Z 2 \ Zx),H(V 2 \ Vi)} - H{V 2 \ Y, V±) 



= H{Z 2 | Zi) - (V2 | Y, Vi) = H(V|17) - H(V|y) 



= /(y ; v|t/) 



(51) 



This is one of the corner points of the rate region given in equation (|50[) . Choosing the permutation tta(') to be 
the derangement gives us the other corner point and time sharing between the two points yields the entire rate 
region of equation (|50[) . The rate needed to reconstruct U,V at the decoder coincides with the Berger-Tung 



We note that this implies that our theorem recovers the rate regions of the problems considered by Wyner 
and Ziv [44], Ahlswede-Korner- Wyner [42,63], Berger and Yeung [47] and Slepian and Wolf [62] since the 
Berger-Tung problem encompasses all these problems as special cases. 

7.4 Lossless Reconstruction of Modulo-2 Sum of Binary Sources 

In this section, we show that Theorem [T] recovers the rate region derived by Korner and Marton [39] for the 
reconstruction of the modulo-2 sum of two binary sources. Let X, Y be correlated binary sources. Let the 
decoder be interested in reconstructing the function F(X,Y) = X (B 2 Y losslessly. In this case, the auxiliary 
random variables can be chosen as U — X, V — Y. Clearly, this function can be embedded in the groups 
Z 2 , Z3, Z4 and 7L 2 ®1 2 . For embedding in 1 2 , the rate region of Theorem [1] reduces to 



It can be verified that embedding in Z3 or Z4 always gives a worse rate than embedding in Z2. Embedding 
in Z2 © Z2 results in the Slepian- Wolf rate region. Combining these rate regions, we see that a sum rate 
of Ri + R 2 = mm(2H(X ©2 Y),H(X,Y)) is achievable using our coding scheme. This recovers the Korner- 
Marton rate region for this problem [39,57]. Moreover, one can also show that this approach can recover the 
Ahlswcdc-Han rate region [52] for this problem, which is an improvement over the Korner-Marton region. 

8 Examples 

In this section, we consider examples of the coding theorem (Theorem [T]). First we consider the problem of 
losslessly reconstructing a function of correlated quaternary sources. We then derive an achievable rate region 
for the case when the decoder is interested in the modulo-2 sum of two binary sources to within a Hamming 
distortion of D. 



rate region [45,46]. 



□ 



i?i > mm(H(X), H{X ffi 2 Y)), R 2 > mm(H(Y),H(X ® 2 Y)) 



(52) 
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8.1 Lossless Encoding of a Quaternary Function 

Consider the following distributed source coding problem. Let (X, Y) be correlated random variables both taking 
values in Z4. Let X, Z be independent random variables taking values in Z4 according to the distributions Px 
and Pz respectively. Define pi = Px(i),qi — Pz{i) for i = 0, . . . , 3. Assume further that the random variable 
Z is non-redundant, i.e., q± + q 3 > 0. Define the random variable Y as Y — X ©4 Z. Suppose X and Y are 
observed by two separate encoders which communicate their quantized observations to a central decoder. The 
decoder is interested in reconstructing the function Z = (X — Y) mod 4 losslessly. 

Since we are interested in lossless reconstruction, we can choose the auxiliary random variables U, V to be 
U = X,V = Y. The function G(U,V) then reduces to F(X,Y) = (X - Y) mod 4. This function can be 
embedded in several groups with order less than or equal to 16. We claim that this function F(X, Y) can be 
embedded in the groups Z4, Z7, Z2 ©Z2 ©Z2 and Z4 ©Z4 with nontrivial performance. For each of these groups, 
we compute the achievable rate as given by Theorem [T] below. For simplicity, we restrict ourselves to the rate 
regions given by equation (|39|) alone. 

Lets consider the group Z4 first. Define the mappings x = S^' '(x) = x for all x G Z4,?/ = Sy^iy) = 
— y for all y £ Z4 and S F 4 (z) — z for all z 6 Z4. With these mappings, it follows from Definition [5] that 
F(X,Y) is embeddable in Z4 with respect to the distribution Pxy ■ From Theorem [I] it follows that an 
achievable rate region using this embedding is given by 

Rx >max{H(Z),2(H(Z)-H([Z] 1 ))} 

= ma,x{h(q ,q 1 ,q2,q3),2(h(qo,q 1 ,q2,q 3 ) - h(q + q 2 ,qi + q 3 ))} (53) 
R 2 > m^{H(Z),2(H(Z) - H([Z\{))} 

= ma,x{h(q , q u q 2 , q 3 ), 2(h(q Q , q 1 ,q 2 , q 3 ) - h(q + q 2 ,qi + 93))} (54) 

giving a sum rate of 

Rz 4 = Ri + Ri > 2ma,x{h(q (h qi,q 2l q 3 ) 1 2(h(q ,q ll q 2 ,q 3 ) - h(q + q 2 ,qi + q 3 ))} (55) 

It can be verified that F(X, Y) can't be embedded in Z5 or Zg. It can be embedded in Z7 with the following 
mappings. Define x = S x Zv \x) = x for all x 6 Z 4 , y = Sy 7 \y) = —y for all y e Z 4 where — y is the additive 
inverse of y in Z 7 and Sf 7) (0) = 0,^ t) (1) = S ( p v) (A) = l,sf 7) {2) = S { p 7} (5) = 2,sf 7 \i) = sf 7 \<6) = 3. 
Let Z = X ©7 Y. From Theorem [TJ it follows that an achievable rate region using this embedding is given by 



Ri > H(Z) = h(q , (1 -po)qs, (l - po - Pi)q2,p 3 qi,poq 3 , (po+pi)q2, (1 -p 3 )qi) (56) 

R 2 > H[Z) = h(qo, (1 -po)qs, (1 - Po - Pi)q2,p 3 qi,Poq 3 , (po+Pi)q2, (1 -^3)51) (57) 

giving a sum rate of 

Rz 7 = Ri +R 2 > 2h(q , (1 -po)q 3 , (1 -p - Pi)q2,p 3 qi,Poq 3 , {po+Pi)q2, (1 ~ p 3 )qi) (58) 



Of the three abelian groups of order 8, it can be verified that embedding F(X, Y) in Zg results in the same 
rate region as given by equations ([56]) and (|57[) and embedding F(X, Y) in Z2 © Z4 results in the same rate 
region as given by equations (|53|) and (f54|) . So, we consider embedding F(X,Y) in Z2 © Z2 © Z2. Recall 
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(a) Sx(-) 



(b) Sy(-) 



(c) S F (-) 



X 


X 





000 


1 


001 


2 


100 


3 


101 



Y 


Y 





000 


1 


010 


2 


100 


3 


110 



z 


S F (z) 


000,011 





001,110 


1 


100,111 


2 


010,101 


3 



Table 1: Mappings for embedding F(X, Y) in Z 2 © Z 2 © Z 2 



that elements of the abelian group Z 2 © Z 2 © Z 2 can be treated as 3 bit vectors over Z 2 . The mappings 

s ^m.m.) { ^ s (z 2 e Z2 eE 2)(0 and s iz M ^)^ are as given in Tableffl 

Define the random variable Z = U ffi V where © is addition in Z 2 © Z 2 © Z 2 . With these mappings, an 
achievable rate region can be derived using Theorem [T] as below. Choose the permutation 7i"z 2 ©z 2 ©z 2 (■) as 
7r(l) = 2, 7r(2) = 3, 7r(3) = 1. Encoding is carried out in 3 stages with the corresponding rates being 

Rn = 0,R 21 = H(Z 2 ) (59) 
#12 = H(Z 3 | Z 2 ),R 22 = (60) 
R 13 = H{Z X | Z 2 ,Z 3 ), R 23 = R{Z X | Z 2 , Z 3 ). (61) 

Summing over the 3 stages of encoding, we get an achievable sum rate of Ri+R 2 > H(Z) + H(Z\ \ Z 2 , Z 3 ) = 
2H(Z) — H(Z 2 . Z 3 ). In terms of pi, qi, this sum rate can be expressed as 

-Rz 2 ©z 2 ©z2 - R i + R 2 > 2Wpo 2 qo,Pi3q3,Po 2 qi,Pi3qo,Po 2 q 2 ,Pi3qi,Po 2 q3,Pi3q 2 ) - Hpo 2 qo 2 ,Pi3qi3,Po 2 qi3,Pi3qo 2 ) 

(62) 

where p 02 = p Q + p 2 ,pi 3 = Pi + P3, <702 = qo + q 2 and qi 3 = qi + q 3 . 

Embedding F(X, Y) in groups of order 9 to 15 result in rate regions which are worse than the ones already 
derived. We next present an achievable rate region when F(X, Y) is embedded in Z4 ©Z4. We use the mappings 
Sf 40Z4) (x) = xO for all x £ Z 4 , S^®^ \y) = Oy for all y £ Z 4 and Sf (xy) = (x,y) for all {x,y) G 1%. 
This embedding corresponds to reconstructing the sources X and Y losslessly and the rate region coincides with 
the Slepian-Wolf rate region. 

Rz 4 ®z, = Ri + R 2 > H{X,Y) = H(X) + H(Z) = h(p ,Pi,p 2 ,p 3 ) + h(q ,qi,q 2 ,q 3 ) (63) 

Combining equations (f55|) . ([58]) . ([62)1 and (|63|) gives us an achievable rate region for this problem. Each of 
these achievable rate regions outperform the others for certain values of Px and Pz ■ This is illustrated in Table 

m 



8.2 Lossy Reconstruction of the Modulo-2 Sum of Binary Sources 

This example concerns the reconstruction of the binary XOR function with the Hamming distortion criterion. 
The rate region of Theorem [T] is very cumbersome to calculate analytically in the general case. So, we restrict 
our attention to the case of symmetric source distribution and additive test channels in the derivation below 
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Px 


P Z 


Pz 4 


Rz 7 




Pz 4 ©Z 4 


U 4 4 4 J 


ri i il 
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rl 1 1 1 i 
I3 10 2 15 J 


r3 1 1 2i 
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3.6847 


4.5925 
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3.4633 


\— — — — 1 

UO 30 30 30 J 


[331 1 1 
120 4 20 20 J 


2.308 


2.7065 


1.9395 


1.7815 



Table 2: Example distributions for which embedding in a given group gives the lowest sum rate. 

where the intention is to demonstrate the analytical evaluation of the rate region of Theorem[T] We then present 
plots where the entire sum rate-distortion region is computed without any restrictive assumptions. 

Consider a binary correlated source (X, Y) with symmetric joint distribution Pxy(0, 0) = Pxy (1, 1) = s/2 
and Pxr(l,0) = Pxy (0,1) = p/2. Suppose we are interested in reconstructing F(X,Y) = X ffi 2 Y within 
Hamming distortion D. We present an achievable rate pair for this problem based on Theorem [1] and compare 
it to the achievable rate region presented in Lemma [T] It was shown in [64] that it suffices to restrict the 
cardinalities of the auxiliary random variables U and V to the cardinalities of their respective source alphabets 
in order to compute the Berger-Tung rate region. Since the scheme presented in Lemma [T] is based on the 
Bcrger-Tung coding scheme, the rate region 1ZT>bt for this problem can be computed by using binary auxiliary 
random variables. 

Let us now evaluate the rate region provided by Theorem [1] for this problem. The auxiliary random variables 
U and V are binary and suppose the test channel PxyPu\xPv\y is fixed. The function G(U,V) which is the 
optimal reconstruction of X ©2 Y given U and V can then be computed. In general, this function can take any 
of the 16 possible values depending upon the test channel PxyP\j\xPv\y ■ 

Let us choose the auxiliary random variables U and V to be binary and for ease of exposition, let them 
be defined as U — X ©2 Qi and V = Y ©2 Q2- Here Qi,Q2 are independent binary random variables with 
P(Qi = 0) = qi,i = 1,2. Let pi = 1 — qi,i — 1,2. Define a — q\q2 + P1P2 and j3 = 1 — a. Once the test 
channel PxyPjj\xPv\y 1S thus fixed, the optimal reconstruction function G(U, V) that minimizes the probability 
P(F(X, Y) ^ G(U, V)) can be computed. It can be showed that 



G(U,V) 





U@ 2 V 

u® 2 v 
1 



a > p, a < q 
a > p,a > q 
a < p,a < q 
a < p,a > q 



(64) 



where a denotes the complement of the bit a. The corresponding distortion for these reconstructions can be 
calculated as 



D(a) = < 



P 
& 
a 

q 



a > p, a < q 
a > p, a > q 
a < p, a < q 
a < p, a > q 



(65) 



Clearly, no rate need be expended if the function to be reconstructed is G(U, V) = or G(U, V) = 1. It is 
also easy to see that the rates needed would be the same for both G(U, V) — U (B2V and G(U, V) — U ©2 V. Let 
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us therefore consider only reconstructing G(U, V) = U® 2 V . It can be shown that this function is embeddable in 
■oup, 



the groups Z 2 , Z 3 , Z4 and Z 2 ©Z 2 . Let us consider the group A = Z 2 . The associated mappings S{f\-), Sy^(-) 



and Sq (•) are all identity mappings. In this case, we have only one digit to encode. Further, note that 
P{Z X = 0) = P{Ux © 2 Vi = 0) = qa + 

The rates of the encoders are given by equations ([55)) and ([^0]) to be 



flu = min{JI(?7 1 ), ff(Zi)} - ff(E^i | X) 

= min{l, /i(<7a + p/3)} — h(qi) 

= h(qa+pf3) - h(qi) 
R21 = min{fl-(Ki),ff(Zi)} - fr(Vi I K) 

= min{l, h(qa + p/3)} - h(q 2 ) 

= h(qa+pf3) - h(q 2 ) 

where h(-) is the binary entropy function. Thus, an achievable rate region for this problem is 



(66) 



(67) 



K= (J {(Ri,R2,D): Ri > h(qa+ P f3) -h(qi),R 2 > h{qa + p/3) - % 2 ), D > D(a)} (68) 

0<9l,92<l 

where D(a) is given in equation (p5|l . Rate points achieved by embedding the function in the abelian groups 
Z3 , Z4 are strictly worse than that achieved by embedding the function in Z 2 while embedding in Z 2 © Z 2 gives 
the Slepian-Wolf rate region for the lossless reconstruction of (X, Y) . 

We now plot the entire sum rate-distortion region for the case of a general source distribution and general 
test channels Pjj\Xi Pv\y an d compare it with the Berger-Tung rate region 1ZT>bt of Fact[T] 

Comparison of sum rate-distortion regions of the two coding schemes 




0.15 0.2 

Distortion D 



Figure 1: Sum rate-distortion region for the distribution given in Table [3] 
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Lower convex envelope of the sum rate-distortion regions of the two coding schemes 
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Figure 2: Lower Convex envelope of the sum rate-distortion region 



PXY 





1 



1 


0.3381 
0.2291 


0.1494 
0.2834 



Table 3: Joint distribution used for example in Figures [Tl2l 

Figures [1] and [2] demonstrate that the sum rate-distortion regions of Theorem [1] and FactfTJ Theorem [1] offers 
improvements over the rate region of Fact [1] for low distortions as shown more clearly in Figure [3] The joint 
distribution of the sources used in this example is given in Table [3l 

Motivation of choosing this example is as follows. Evaluation of the Berger-tung rate region is a computa- 
tionally intensive operation since it involves solving a nonconvex optimization problem. The only procedure that 
we are aware of for this is using linear programming followed by quantizing the probability space and searching 
for optimum values [64] . The computational complexity increases dramatically as the size of the alphabet of 
the sources goes beyond two. Hence we chose the simplest nontrivial lossy example to make the point. But 
then we have to deal with the fact that there are only 3 abelian groups of order less than or equal to 4. One 
of the three groups corresponds to the Berger-Tung bound. We would like to remark that even for this simple 
example, the Berger-Tung bound is not tight. We expect the gains afforded by Theorem [T] over the rate region 
of Lemma Q] would increase as we increase the cardinality of the source alphabets. 
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Comparison of the two coding schemes for low distortion 



Comparison of the two lower convex envelopes for low distortion 
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Figure 3: Zoomed versions of Figures [T] and [2] 

9 Conclusion 

We have introduced structured codes built over arbitrary abelian groups for lossless and lossy source coding 
and derived their performance limits. We also derived a coding theorem based on nested group codes for recon- 
structing an arbitrary function of the sources based on a fidelity criterion. The encoding proceeds sequentially in 
stages based on the primary cyclic decomposition of the underlying abelian group. This coding scheme recovers 
the known rate regions of many distributed source coding problems while presenting new rate regions to others. 
The usefulness of the scheme is demonstrated with both lossless and lossy examples. 



Acknowledgements 

The authors would like to thank Professor Hans-Andrea Locligcr of ETH, Zurich and Dr. Soumya Jana of 
University of Illinois, Urbana-Champaign for helpful discussions. 



Appendix 

A Good Group Channel Codes 

We prove the existence of channel codes built over the space Z^ r which are good for the triple (Z,S,Pzs) 
according to Definition [lOl Recall that the group Z p >- has (r — 1) non-trivial subgroups, namely p J Z p >-, 1 < i < 
r — 1. Let the random variable Z take values from the group Z pr , i.e., Z = Z pI . and further let it be non- 
redundant. Let Hom(Z™ r ,Zp r ) be the set of all homomorphisms from Z™ r to Zpv. Let </>(■) be a homomorphism 
picked at random with uniform probability from Hom(Z™ r , Zpv). 
We start by proving a couple of lemmas. 
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Lemma 4. For a homomorphism <p(-) randomly chosen from Hom(Zp r , Zpv ) , the probability that a given se- 
quence z n belongs to ker(<^) in Z™ r depends on which subgroup of Z pr the sequence z n belongs to. Specifically 

PMz") = 0*) = { P ~ ir ~ i)k lfzn 6 ^A^,0 < i < r 

[ l !jz»e/z; r 

Proof: Clearly, z n G p r Z p V implies z™ = 0™[f|. In this case, the probability of the event {(p{z n ) = 0} is 1. 

Let the k x n matrix <£> be the matrix representation of the homomorphism </>(•). Let the first row of $ 
be (ai, . . . , a n ). Consider 4>i : Z pr — » Z pr , the homomorphism corresponding to the first row of $. The total 
number of possibilities for </>i(0 is {p r ) n . 

Let us consider the case where z n G Z p r\pZ p ,.. In this case, z n contains at least one element, say Zi which 
is invertible in Z p r. Let us count the number of homomorphisms </>(•) that map such a sequence z n to a given 
c G Zpv. We need to choose the k homomorphisms </>i(-),l < i < k such that 4>i(z n ) — Ci for 1 < i < k. Let 
us count the number of homomorphisms ^i(-) that map z n to c\. In this case, we can choose ctj,j i to be 
arbitrary and fix on as 

V j?i / 

Thus the number of favorable homomorphisms </>i(-) is (p r )( n_1 ). Thus the probability that a randomly chosen 
homomorphism <j)\{-) maps z n to c\ is Since each of the k homomorphisms 4>i can be chosen independently, 
we have 

p{4>{z n ) = c) = P ~ rk if 2 " g z;wz p v (71) 

Putting c = fc in equation (jTT]) , we see that the claim in Lemma 2] is valid for z n G Z pr \pZ pr . Now, consider 
z™ G p 4 Z™ r \p l+1 Zp r for a general < i < r. Any such z n can be written as p % z n for z n G Z pr .\pZ p ,.. Thus, the 
event {4>(z n ) = 0} will be true if and only if {<f>(z n ) = t} for some t G p r ~ l Z k )r . Hence, 



/ 



P((f>(z n ) = 0) = P 



U wn = t)\ (72) 

= P ^ n )=t) (73) 

f p r 

= \p r - l 1 k pr \p- rk (74) 

= p-^ k (75) 

This proves the claim of Lemma [4] □ 
We now estimate the size of the intersection of the conditionally typical set A™(s n ) with cosets of p l Zp r in 



3 If we consider homomorphisms from Z™,- to for an arbitrary integer m, all such homomorphisms have dZ™,. as their kernel 
where d = (p r ,m) is the greatest common divisor of p T and m. Unless d = p r , there would be exponentially many z n for which 



2* 

and all such m give identical performances as m = p 



P{4>(z n ) = 0) = 1 for all <j> S Hom(Z™ r , Zjjj) and this results in bad channel codes (see equation (11071 ). Thus, p r has to divide m 
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Lemma 5. For a given z n G A™(s n ), consider (z n + p l Z p V) , the coset of p ir L pr inZ pr . Define the set Si^(z n , s n ) 
as Si. e (z n , s n ) = [z n + p ir L pr ) n A"(s"). A uniform bound on the cardinality of this set is given by 

^]Qg\Si, e \<H(Z\S)-H([Z\i\S)+5(e) forO<i<r (76) 

where 5(e) — > as e — > 0. T/ie random variable [Z]i is defined in the following manner: It takes values from the 
set of all distinct cosets of p l Z p r in Z p r. The probability that [Z]» takes a particular coset as its value is equal 
to the sum of the probabilities of the elements forming that coset. 

P([Z]i = a + p*Z pr | S = s)= P zls (z | s) V S G S. (77) 

zea+p if l p r 

We have the nesting relation S i+ i^(z n ,s n ) C Si_ e (z n , s") for < i < r — 1. However, each nested set is 
exponentially smaller in size since H([Z]i) increases monotonically with i. Thus, with the same definitions as 
above, we also have that 

l - log (\S ite (z n , s n )\ - \S i+he (z n , s n )\) < H(Z\S) - H([Z]i\S) + Si(e) for0<i<r-l (78) 
where <5i(e) ^0 as e — > 0. 

Proof: The set Si i€ (z n ,s n ) can be thought of as all those sequences z n € A"(s") such that the difference 
w n = z™ — z™ e p l ^ p r. Let be a random variable taking values in j/Z pr - and jointly distributed with (Z, S 1 ) 
according to P\y\zs- Dehne the random variable Z = Z + W. Let Py/\zs be such that P z ~ s = Pzs- Then, for 
a given distribution Pw\zs> every sequence z n that belongs to the set of conditionally typical sequences given 
(z n ,s n ) will belong to the set Si te (z n ,s n ). Conversely, following the type counting lemma and the continuity 
of entropy as a function of probability distributions [57], every sequence z n <G S^ e (z n ,s n ) belongs to the set 
of conditionally typical sequences given (z n ,s n ) for some such joint distribution Pw\zs- Thus estimating the 
size of the set Si^ e (z n , s n ) reduces to estimating the maximum of H(Z \ Z, S), or equivalently the maximum of 
H (Z, W | S) over all joint distributions Pzsw such that P(z+w),s = Pzs- 

We formulate this problem as a convex optimization problem in the following manner. Recall that the 
alphabet of Z is the group Z pr . Hence, H(Z,W \ S) is a concave function of the \Z\ \S\ \p l 7h p r \ variables 
Pzsw{Z — z,S = s,W = w),z G Z,S G S, w G p % 7h p r and maximizing this conditional entropy is a convex 
minimization problem. Since the distribution Pzs is fixed, these variables satisfy the marginal constraint 

Pzw\s(Z = z,W = w\S = S ) = P zls (z\s) VzgZ,,sg5 (79) 

The other constraint to be satisfied is that the random variable Z = Z + W is jointly distributed with S in the 
same way as Z, i.e., P zs = Pzs- This can be expressed as 

Pzw\s(Z = z-w,W = w\S = s) = P z]s (z | s) VzeZ,seS. (80) 
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Thus the convex optimization problem can be stated as 
minimize — H(Z, W \ S) 

subject to ^2 Pzw\s(z,w | s) = P z \s(z s) Vz G Z, s G S, 

E p zw\s{z - w,w I s) = P z \ s {z | s) VzeZ,seS. (81) 

Note that the objective function to be minimized is convex and the constraints of equations (|79[) and ([50)1 on 
P ZW \ S (Z — z, W = w | S = s) are affine. Thus, the Karush-Kuhn- Tucker (KKT) conditions [65] are necessary 
and sufficient for the points to be primal and dual optimal. We now derive the KKT conditions for this problem. 
We formulate the dual problem as 

D{Pzw\s) = ~ E P sC<0 E p zw\s(z, w | a) log \ 

tts \zez,^ P%r P zw \s(z,w\s) 

+ 2\I Xz s \ E p zw\s(z -w,w\s)- P z \s(z I s) 

+ E ^ ( E P zw\s(z,w\s)-P z \s(z\s)\ (82) 
zez.ses y u ;ep i z pr / 

where {X z . s }, {lz,s} are the Lagrange multipliers. Differentiating with respect to Pzw\s{Z = z,W = w \ S — s) 
and setting the derivative to 0, we get 

dD(P Z w\s) = ps(s)(1 +logP {ZjW | s)) + A + = o (83) 

<Jr Z w\S\ z -> w I s ) 

=> \z+w),s + jz,s = -Ps(s)(l + logP Z w\s(z,w | s)) Vz e Z,s e 5, to e p l Z pr . (84) 

Summing over all z G p'Zy. for a given s E S, we see that for all to € _p l Z pI . , the summation 5Z z6p » z [\z+w),s + 7z,s) 
is the same. This implies that 

TQ Pzw\s( z i w I s ) = constant Vto G p l ^ p ,-, Vs G 5. (85) 

These |6>|p r_4 equations form the KKT equations and any solution that satisfies equations (|79[) . (|50"|) and ([55]) is 
the optimal solution to the optimization problem (|81|) . We claim that the solution to this system of equations 
is given by 

p i i ^ Pzjsjz I s)Pz|s(z + w I s) 

PZ\S[Z + P'Zpr I S) 

For this choice of Pzw\si z i w I s )> we now show that equation (|T9[) is satisfied. 

V- P / I \ p z\s(z I s)P Z | S (^ + w | s) 

~i rri Pz\s{z + P l V s) 



2^ Pz|s(z + w|s) 



P z \ s (z\s) VzgZ,sg5. (89) 
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Next, lets show that the choice of P zw \g(z,w \ s) in equation (|86p satisfies equation (|80[) . 

p 17 xv I q \ P z \ s (z-w | | a) 

e ? a ^Wi^ = *-^ = «'|5 = -)= E p^^TT^TR (90) 

p f \ \ P z \s(z-w | a) 

= Pz\s{z\s) VzeZ,seS. (92) 
Finally, we show that this choice of Pzw\s{ z , w I s ) satisfies the KKT conditions given by equation (|85|) . 



P zw \ s (z,w\s) = II — 1 — 93 



| s)P Z |g(z + w | s) 
P Z | S (^+P 1 V I a) 

, \ i r V TT ^Iis(^ I «) (94) 

p Z i S (p 1 V I s )/ ii 1 



which is independent of w and is the same for any w 6 p l Z pr . Thus, equation (|86p indeed is the solution to the 
optimization problem described by equation (|8ip . Let us now compute the maximum value that the entropy 
H(W | Z, S) takes for this choice of the conditional distribution P Z \y\S- 

TTsltz \ w ^z pr P wlzs (w\z,s)J 

„ , J Pz\s(z + w\s) P zls (z+p l Z pr | s) \ 

Let PC be the set of all distinct cosets of p*Z p ,. in Z pr and let T>C(z) be the unique set in VC that contains z. 
Let us evaluate the summation in the brackets of equation (|96p first. 

v P Z |s(^ + ^|s) P Z | 5 (z+p*Z p ,. | g) ^ P Z | S (z + W | g) P Z | S (PC(z) | g) 

P zls (z + p% r |a) ° S P Z | S (z + W |s) 2jj P z , s (PC(z) I «) ° g P z|s (^ + w;| s ) 1 j 

= logP z , s (PC(z) | a) + L Pz|s(z |s) 

P Z \s{T>C(z) | s) 

(98) 

This sum is dependent on z only through the coset T>C{z) to which z belongs. Thus, the sum is the same for 
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any two z that belong to the same coset of p l Z pr in Z pr . Thus, we have 



H(W | Z,S) = ^Ps(s) E E P ^ 2 I s ) l0 S^|5(T | s) 



seS Tevc zeT 



Pz\s(T | a) 



E E logP Z | S (T 



seSTevc 



E z >eT P z\s(z' | g )log Pz|s ^ |s) 

fWr I ^) 



- (Z | S) - H([Z] t | S) 
where \Z\ is as defined in Lemma [5l 



s£S 



Pz\s(T | a) 



(99) 

(100) 

(101) 

(102) 
□ 



We are now ready to prove the existence of good group channel codes. Let Z take values in the group 
Z pr and further be non-redundant. Coding is done in blocks of length n. We show the existence of a good 
channel code by averaging the probability of a decoding error over all possible choices of the homomorphism 
</)(■) from the family Hom(Z" r ,Zpr). Let H be the parity check matrix and C be the kernel of a randomly chosen 
homomorphism </>(•). 

The probability of the set B e (C) can be written as 



P zs (B e (C))= E Pzs(z n ,s n )I 



\ 



(8»,«")gAJ(Z,S) 



(103) 



/ 



where /(-E) is the indicator of the event E. Taking the expectation of this probability, we get 



E(P ZS (B £ (C))) = E Pzs(z n ,s n )P 



|J (cf>(z n ) = t(z n )) 



(S»,s»)6A»(Z,S) 



/ 



< E Pzs(z n ,s n )P 

(z n ,s n )eA^(Z,S) 



< E Pzs(z n ,s n )P 



|J = ct>(z n )) 



J 



( 



\ 



|J {cj>{z n - z n ) = fe ) 



(f", s ™)eA"(z,s) 



E 

(z'V")£A"(Z,S) 



(104) 

Pzs(z n ,s n ) 
(105) 
(106) 



/ 



where 5% — > as n — ► oo. 

We now derive a uniform bound for the probability that for a given (z™, s n ) £ A™(Z, S), a randomly chosen 
homomorphism maps z n to the same syndrome as z 11 for some z n such that (z n , s n ) <E A™(Z, S). From Lemma 
|4] and El w e see that this probability depends on which of the sets Si te (z n , s n ), < i < r the sequence z n belongs 
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to. 



u 



{z n - z n ) = k 



. k i«, s ")eAj(z,s) 



< 



E 



P{(j>{z n - z n ) = k ) 



= E E p (^ r 

i=0 S~eS i ,«(«»,a»)\Si+i,«(a",s n ) 
(a)^ |5,, e (z n , S ")| - |5 l+1 . e (z", S ")| 



z") = fe ) 



E 

< 2^ ex P2 

i=0 



p(r-i)k 



H{Z\S) - H([Z\i\S) - -(r - i) logp + fc(e) 



(107) 

(108) 
(109) 
(110) 



where (a) follows from Lemma|4]and (b) follows from LemmaO If this summation were to go to zero with block 
length, it would follow from equation (|106|) that the expected probability of the set B e (C) also goes to zero. 
This implies the existence of at least one homomorphism </>(•) such that the associated codebook C satisfies for 
a given e > 0, Pzs(B € (C)) < e for sufficiently large block length. 

The summation in equation (|110[) goes to zero if each of the terms goes to zero. This happens if 



k r 



or equivalently 



fc(» 



logp r >H(Z\S)-H([Z]i\S)+S 2 (e) for % = 0, . 



■ log p r > max 



(H(Z\S)-H([Z] i \S)) + S 2 (e) 



(111) 



(112) 



n ' o<i<r \r — i 

It is clear that in the limit as n — *■ oo, good group channel codes exist such that the dimensions of the associated 
parity check matrices satisfy equation (j2"?| . When C is a good channel code, define the decoding function 
i> : Z£ r x S n -> %%r for a given (z n , s n ) as the unique member of the set {z n : Hz n = Hz n , (z n , s n ) G A£(Z, S)}. 



B Good Group Source Codes 

We prove the existence of source codes built over the space Z™ r which are good for the triple (X,U,Pxu) 
according to Definition [9] Let the random variable U take values from the group Z pr , i.e., U — Z pr and let U 
be non-redundant. Let <fi: Z™ r — > Z k r be a homomorphism for some k to be fixed later. The codebook C is the 
kernel ker(^) of this homomorphism. Note that ker(</>) < Z™ r and hence the codebook has a group structure. 
We show the existence of a good code C by averaging the probability of error over all possible choices of </>(•) 
from the family of all homomorphisms Hom(Z™ r , Z* r ) . 

Recall the definition of the set A e {C) from equation (|22|) . The probability of this set can be written as 

p(mc)) = Px{* n )i [ U (* n > u ") e u n ( 113 ) 

at™ \u n eC ) 
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The expected value of this probability is 

(114) 



E(P(A e (C))) = £Px(*")P ( |J (x n ,u n ) £ A-(X,U)) 



(115) 



For a typical x™, let us compute the probability that there exists no u n € C jointly typical with the source 
sequence x n . Define the random variable 9 (x n ) as 

Q(x n )= ]T Uu^ey. (116) 

<d(x n ) counts the number of u n sequences in the codebook C that are jointly typical with x n . The error event 
E given that the source sequence is x n is equivalent to the event {<d(x n ) = 0}. Thus, we need to evaluate 
the probability of this event. Note that Q(x n ) is a sum of indicator random variables some of which might be 
dependent. This dependence arises from the structural constraint on the codebook C. For example, u" G C 
implies that ku™ € C as well for any k S Z pr . We use Suen's inequality [61] to bound this probability. 

In order to use Suen's inequality, we need to form the dependency graph between these indicator random 
variables. We do this in a series of lemmas. We first evaluate the probability that a given typical sequence 
belongs to the kernel of a randomly chosen homomorphism. Since U is assumed to be non-redundant, by Lemma 
[H we have 

P(u n eC) = p - rk (117) 

We now turn our attention to pairwise relations between the indicator random variables. For two n-length 
sequences u™, u 2 , define the matrices Mfc^u™, u 2 ), 1 < k,l < n and k ^ I as 



M M «, M ") 



U 2 k U 2 l 



(118) 



Let m,k,i(ui,u 2 ) be the determinant of the matrix Mkjiui ,u 2 ). Define the set 

MK,<) 4 {m M K, t#) : u" 1 exists} (119) 

Note that the set M(ui,u 2 ) is non-empty since is assumed to be a non-redundant sequence. Let D(u",u 2 ) 
be the smallest subgroup of Z pr that contains the set M(u", u 2 ). As will be shown, the probability that both u" 
and u 2 belong to the kernel of a randomly chosen homomorphism depends on -D(u™, u 2 )■ For ease of notation, 
we suppress the dependence of the various quantities on the sequences u™,u 2 m what follows. 

Lemma 6. For two non-redundant sequences u™, u 2 , the probability that a random homomorphism 4>: Z™ r — ► 
maps the sequences to k is 

P(<£«) - 0K) - fe ) = p -(*-«)fc if £)«, v») = p% r , < i < r (120) 

Proof: Let the homomorphism </)(•) be decomposed as <^>i : Z™ r — » Z pr , 1 < i < k. We first count the number 
of homomorphisms <fii(-) that map both u" and u 2 to 0. Recall that </>i(m") can be expressed as the linear 
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combination = Xw=i a j u ij f° r a j ^ ^V'j 1 ^ J — n - Thus, we need to find the number of solutions 

that simultaneously satisfy the equations 

n 
n 

^a iW2i =0 (122) 

If D(ui,U2) — p'^pr, then there exists some 1 < k < n such that u^. exists and m^^* € p I Z pr \p i+1 Z pI . 
for some 1 < j* < n, _?'* 7^ fc. Fix such a /c. Then, any solution to the equation (|12ip must be of the form 
otj,j 7^ k arbitrary, otu = Sj^fe a j u ij f° r some k such that u^ 1 exists. Thus, the total number of solutions 

to equation (|12ip is p r ( n ~ x \ Substituting one such solution into equation (|122p , we get 

<XjU 2 j = a 3 u V ~ u ik u ^ Yl a i Ul i ( 123 ) 



^ Otj (UlkU2j ~U2kUlj) (124) 



Of the p r (™ -1 ) choices for {ai}™ =1 , we need to find those that satisfy Sj^fc a j m k,j — 0. We allow aj to be 
arbitrary for j 7^ k,j* and solve the equation aj-m^.j* = — X}jyfc i» a j m k,j- It is clear that the summation in 
the right hand side yields a sum that belongs to p l 7L pT . Since k, j* are chosen such that nikj* G p i Z pr \p l+1 l lpr , 
it follows from Lemma|S]in Appendix iDl that this equation has p l solutions for Oj« for each of the p r ( n_2 ) choices 
of aj,j 7^ k,j*. Once etj,j 7^ k is fixed, afc is automatically fixed at au = — u^ k X^j^fc a j u ij- Thus, the total 
number of solutions that simultaneously satisfy equations (I12ip and (|122p is p l p r ^ n ~ 2 \ 

It follows that the probability of a randomly chosen homomorphism </>i(-) mapping both u™, to is given 
by p l /p 2r . Since each of the k homomorphisms 1 < i < k can be chosen independently, we have 

P(0K) = 0K) = 0) = p-^)* ( i26) 
when D(u™, uV,) — p % 7L p r for some < i < r. This proves the claim of Lemma |5] □ 

Suppose Ui and are non-redundant sequences. It follows from Lemmas |4] and [6] that the events l{ u ™gc} 
and l{ u n 6 cj are independent when D(u™,u%) = Z pr . In order to infer the dependency graph of the indicator 
random variables in equation (|116p . we need to count the number of sequences v% for a given w" such that 
D(u™, u 2 l ) — p l ^ p r for a given 1 < i < r. This is the content of the next lemma. 

Lemma 7. Let it™ be a non-redundant sequence. Let Di(tt™), < i < r be the set of all u% sequences such that 
D(u r i,U2) = p l 7L pr . The size of 'the set Z?i(u™) is given by 

( rf (V r - i )( n ^ 1 ) — nir-i- 1 )^- 1 )) n < i < r 
|A«)|= r , (127) 

p — 1 i = r 
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Proof: We start by estimating the size of D r (u™), i.e., the set of u% sequences such that D(u^, U2) = 0. Since 
D(u™, u%) = 0, u 2 l must be such that there exists 1 < k < n such that it^ 1 exists and mk.j — for all j =/= k. 
This implies that u\kUij — U2kU\j for all j ^ k. Define rj = u^U2k- It then follows that u-ij — rjuij for all 
1 < j < n which implies that = r\u\ for some r\ 6 Z pr . Since it is assumed that ^ u™, there are p r — 1 
distinct values of r;. Since the sequence it" is non-redundant, it follows that each value of rj results in a distinct 
value of Thus, |D r (w™)| = p r — 1 as claimed in the Lemma. 

Consider the case when D^u^^u^) = V^pr for some < i < r. We count the number of 11% for a given u™ 
such that p l Z pr is the smallest subgroup containing all the set M(u", u^)- Since D(u™, uj) = p l Z p r, U2 must be 
such that there exists 1 < k < n such that exists and rrikj* £ P l ^ P 'Ap t+1 Z pr for some 1 < j* < n, j* ^ fc. 
Consider the matrices Mfc^it™, ), 1 < I < n,l ^ k. Let A^j e < I < n,l =fc k. Fixing the values 

of the determinants mk,i{ui, U2) to be A^, we can solve for the entire sequence wj. Thus, ZA;(ii™) contains 
the union over all permissible values of {A^}^ of those sequences itj such that TOfc^it™, wl?) = Ajy for all 
1 < I < n,l ^ k. 

For a given {&.k,i\i=£k, let us investigate the number of u r 2 l sequences such that m^^u™, u 1 ^ ) = ; for all 
1 < I < n,l =/= k. Consider first the equation mk,i* = &k,i* for some I* ^ k. Since uu is invertible, there are p r 
possible solutions in (u2k,U2i*) for this equation. Now consider the equations rrikj, 1 < I < n, I ^ k, I*. Since 
U2k is already fixed and u\k is invertible, there is precisely one solution to U21 in these equations. Solving these 
(n — 1) equations fixes the sequence u r 2 l . Thus, the number of solutions to for a given it™ and {^k,i\i=jtk is p r ■ 
The number of A k ,i such that {A k ,i}i^ k E p'Z^r 1 is clearly pO"-*)(n-i). For D(v%, v%) = p l Z pr , there must exist 
at least one A fe , ; e p l 1 p r \p l+1 Z pr . The total number of such {A k ,i}i^ k is clearly pO-'K™-!) _ p (r-i-i)(n-i)_ 
Putting these arguments together, we get that the size of .Dj(u™) is p r (p( r_t )( n_1 ) — p( r ~ l ~ 1 )( n - 1 )y This proves 
the claim of Lemma [7] □ 

We are now ready to infer the dependency graph of the indicator random variables in equation (|116j) . The 
number of nodes in the dependency graph is |-A"(x™)|. Let li be the indicator of the event {m™ £ C} and let 
Ii correspond to the ith vertex of the graph. From Lemma [5J it follows that vertices i and j are connected 
(denoted by i ~ j) if ^ Using Lemma [3 the degree of the ith vertex can be bounded by 

p rn — \Dq(ui )| — 1 = j3'+( r_1 )(™ _1 ) — 1. Note that this is an upper bound since not all U2 sequences counted in 
Lemma [7] need belong to A™(x n ). 

One version of Suen's inequality can be stated as follows. Let Ii £ Be(pi),i £ X be a family of Bernoulli 
random variables having a dependency graph L with vertex set X and edge set E(L). Let X = Ei-^i an d 



A = E(X) = J2 lPl . Write * ~ j if £ E(L) and let A = \ E,-; W,-) and * = max < E fc ~<Pfc- Then 



Let us estimate the quantities A, A and 8 for our problem. It follows from equation (|117[) that A = E(0(x")) = 
\A™ (x n )\p~ rk . Uniform upper and lower bounds [57] exist for the size of the set A™(x n ). An upper bound to A 




(128) 
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can be established via Lemmas |5] and [7] as below. 

A = ^EE E ( 7 ^) ( 129 ) 

4EE p WO = ^K) = o) (130) 



= \ E E E p(0«) = 0K) = o) (i3i) 

uj'e^ji") m=i « ;j eA™(x'»)nz3 m («5') 

= 2 E El^OOnAnK)! ^ (132) 

«™GA"(x") m=l ^ ' 



E f(p r -D(i) fc + Ep r t 
i|^(x n )i ( ( P r - 1) E p r (p (r - m)(n - 1] - P ^- m -^ n -^) (Cm ( 134 ) 



< ^ n E ((p- - 1) + E/ (^ ) ! «!•>->) 

/1 ?C a ' ri ) 



where (a) follows from Lemma [6] and (6) follows from Lemma [7] This expression can be further simplified by 
noting that f(m) = \D m (u") decreasing function of m. Thus, the summation in the parentheses 

of equation (fTMj) can be upper bounded by (r - 1)/(1) = (r - l)| J D 1 (u?)|p- fe ( 2r - 1 ). Thus, 

A < \\A^{x n )\ (p r ~ rk + (r- l) P " rW1 - n - 2rfc (l - pM) (135) 

< -|^(a;")|j) r - rfe (l + (r - ljpfr- 1 )^-*- 1 )) (136) 
We now bound the quantity S. 

S = max^E(ij) (137) 

r 

= „max V Y, ^K) = 0) (138) 

< max f p N-(r-l)(n-l) _ ^ p -rfc ( 139 ) 

u™eA™(:r™) V / 

< p r(n-fc)-(n-l) ( 14 q) 

where (a) follows from equation (| 11 T[) and the fact the Pjj\x is a non-redundant distribution. Using these 
bounds, we can bound the terms involved in equation (|128[) . 

\2 I Ant rr n\\2 n -2rk 

A_ > \ A ^ X }} p (141) 

8A - 4\A™(x n )\p r - rk (1 + (r- i) p (r-i)(n-fc-i)) v ' 

|A"(x")|p- r ( fc+1 ) 

- 4(1 + rp('- 1 )("- fc - 1 )) ( ' 142 ' ) 
> \A"(x n )\ Mr _ 1)+k+1) 
~ 8r 
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where the last inequality holds for sufficiently large n. The third term in the exponent in equation (|128p can be 
bounded as 

_A („( r -i)+i) 
65 ~ Q 1 

Combining equations (|143|) and (|144|1 . we get a bound on the probability of the event {Q(x n ) = 0} as 



(144) 



r , (-)„•'<) =,„ ,. N])< j -min ^"^"^ -^, l A "^ n )l p -(»(r-l)+fc+l) i l^(*")l p -(n(r-l) + l) 



(145) 



As long as each of the terms in the minimizations goes to oo as n — > oo, the probability of not finding a jointly 
typical sequence with x n in the codebook C goes to 0. Let x n € A™(X) be a typical sequence. It is well 
known [57] that for sufficiently large n, the size of the set A™(x n ) is lower bounded as 



|^4™(a;™)| > 2 n ( H( J J][X " > ~' Ll ^ 



where ei(e) — > as e — > 0. Therefore, 

\^{x n )\ k 



>-exp 2 (n H(U\X) 
\A n t {x n )\ p - {n[T - l)+k+1) > ± , „ 

j^llp-Wr-D+l) > l exp2 ( /f //(f ,|., Y) 



rfc 



logp 



H(U\X)- (r-l) + 



logp - ei 



V - 11 



fogp - £i 



(146) 

(147) 
(148) 
(149) 



For the probability in equation (|145[) to decay to 0, we need the exponents of these three terms to be positive. 
Equation (|147[) gives us the condition 

k 

- log p r < H(U\X) (150) 
n 

while equations (I148[) and (|149|) together give us the condition 



< - logp r < r(H(U\X) - logy- 1 ) 
n 

Thus, the dimensionality of the parity check matrix satisfies the asymptotic condition 

k ( n ) i„„„r ;„firfTT\v\ „lzTfTT\v\ l„„„r-l| + \ 



lim 



logp r = min(H(U\X),r\H(U\X) - logjf- 1 ^) 



(151) 



(152) 



n — >oo 77, 

where \x\ + = max(x, 0). Combining these results, we see that provided equation (|152|) is satisfied, P(Q(x n ) = 0) 
goes to double exponentially. We now show that there exists at least one codebook C such that the set A t (C) 
has high probability. We do this by calculating the ensemble average of P(A e (C)) over all codebooks C. It 
follows from equation (|115p that 



E(P x (A e (C))) > J2 Px(x n )P(®(* n ^ 0)) 

i"£A?(X) 



(153) 



> (1 - e 2 ) 1 - exp 



\A?(x n )\ , |A^(a;")| ( ?t(T .-i) +fc+ i) |^"(x")| (n _ 1) 
2 8r 6 



(154) 



where e-i — * as rt — > oo. Thus, as long as equation (|152|) is satisfied, the expected value of Px(A e (C)) can be 
made arbitrarily close to 1. This implies that there exists at least one homomorphism such that its kernel is a 
good source code for the triple (X,U,Pxu)- 
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C Good Nested Group Codes 



We now show the existence of good nested group codes satisfying Lemma [3j As was remarked in Definition [8j 
one way to construct a nested group code is to add rows to the parity check matrix of the fine code to get the 
parity check matrix of the coarse code. Let the random variables X, Y, U, V, S be as given in Lemma [3] Let 
the parity check matrices of the codes C\\,C\2 and C2 be Hu, H12 and H2 respectively. Let their corresponding 
dimensions be fcn x n, k\2 X n and k2 x n respectively. In order to ensure nesting, impose the following structural 
constraints on these matrices. 



H12 = 





, H 2 = 


H\2 






. AH2 . 



(155) 



Let the dimensions kn,ki2 and k 2 satisfy equations P§j) - (p?Tj) . 

Generate random H2,H\2 matrices by constructing the matrices Hu, AH±, AH2 independently by picking 
entries uniformly and independently from the group T, r . From the proofs in Appendices [A] and [B] it follows 
that the codes C\\,C\2 and C2 are with high probability good source and channel codes respectively for the 
appropriate triples. By union bound, it follows then that there exists a choice of Hu, AH\ and AH2 such that 
the codebook C2 is a good channel code and the nested codes Cn and C12 are simultaneously good source codes 
for their respective triples. This proves the existence of good nested group codes as claimed in Lemma [3] 



D Linear Equations in Groups 

We now present a lemma on the number of solutions over the group Z pr for a linear equation in one variable. 

Lemma 8. Let a G p ir L p r\p lJrX TL pT for some < i < r. Then, the linear equation ax — b has a solution in x if 
and only if b £ p l Z pr . In that case, there are p 1 distinct solutions for x over the group 1 pr . 

Proof: It is clear that the equation ax = b cannot have a solution if b ^ p l 7L p r. The rest of the proof proceeds 
in two stages. We first show that if there exists at least one solution to the equation ax = b, then there exists 
p l distinct solutions. We then show that at least one solution exists for every b e p l Z pr . Together, these imply 
Lemma |SJ 

Suppose there exists at least one solution x\ to the equation ax = b. Then, for any t € p r ~ % 'L pr , x\ +t is also 
a solution and all such solutions are distinct. Conversely, if x\, X2 are both solutions, then x\ — X2 £ p r ~ l 1i pr . 
Thus, existence of at least one solution implies the existence of exactly p 1 solutions. Now consider the number 
of distinct values of the set {ax: x £ Z pr }. Since every distinct value repeats itself exactly p l times and there 
are p r elements in this set, it follows that the number of distinct values is p r ~ % . This is exactly the size of the 
subgroup p t Z pr which implies that ax = b has exactly p l solutions for every element b G p % 7L pT . □ 

E T is non-empty 

Recall the definition of T from Section[6]as T = {A : A is abelian, \Q\ < \A\ < a/3, G(U, V) C A with respect to Pjjv}- 
Let \U\ = a, |V| = 0. We now show that the function G(U, V) can always be embedded in some abelian group 
belonging to T. Consider the function G\(U, V) = (U, V). Clearly, Gi(U, V) C Z a © Zp with respect to Puv 
for any distribution Puv- Since there is an obvious surjective mapping between the functions G\{U,V) and 
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G(U, V), it follows from Definition [5] that G(U, V) C Z Q © with respect to Puv- Since |Z a © Z^l = a[3, it 
follows that this group belongs to the set T and hence T is always non-empty. 

References 



[9 
[10 

[11 
[12 

[is; 

[14 
[15 
[16 



P. Elias, "Coding for noisy channels", IRE Conv. Record, part. 4, pp. 37-46, 1955. 

T. J. Goblick, Jr., "Coding for a discrete information source with a distortion measure" , Ph.D. dissertation, 
Dept. Electr. Eng., MIT , Cambridge, MA, 1962 

R. L. Dobrushin, "Asymptotic optimality of group and systematic codes for some channels", Theor. Probab. 
Appl, vol. 8, pp. 52-66, 1963. 

R. G. Gallager, Information Theory and Reliable Communication. John Wiley and Sons, Inc., 1968. 

H. Yamamoto and K. Itoh, "Source coding theory for multiterminal communication systems with a remote 
source," The Transactions of the IECE of Japan, vol. E-63, pp. 700-706, October 1980. 

T. J. Flynn and R. M. Gray, "Encoding of correlated observations," IEEE Trans. Inform. Theory, vol. IT- 
33, pp. 773-787, November 1987. 

A. Orlitsky and J. R. Roche, "Coding for computing," IEEE Trans. Inform. Theory, vol. IT-47, pp. 903-917, 
March 2001. 

N. Ma and P. Ishwar, "Two-terminal distributed source coding with alternating messages for function 
computation", Proc. IEEE Interational Symposium on Inform. Theory, Toronto, Canada, 2008. 

A. Giridhar and P. R. Kumar, "Computing and communication functions over sensor networks" , IEEE 
Journal on selected areas in communications, vol. 23, no. 4, pp. 755-764, April 2005. 

J. Muramatsu and S. Miyake, "Hash property and coding theorems for sparse matrices and maximum- 
likelihood coding", Proc. IEEE International Symposium on Inform. Theory, Toronto, Canada, 2008. 

H. Yamamoto, "Wyner-Ziv theory for a general function of the correlated sources," IEEE Trans. Inform. 
Theory, vol. IT-28, pp. 803-807, September 1982. 

G. Cohen, I. Honkala, S. Lytsyn and A. Lobstein, Covering Codes. North Holland-Elsevier, 1997 

J. H. Conway and N. J. A. Sloane, Sphere Packings, Lattices and Groups. Springer, 1992. 

P. Delsarte and P. M. Piret, "Do most linear codes achieve the Goblick bound on the covering radius?", 
IEEE Trnas. Inform. Theory, vol. IT-32, no. 6, pp. 826-828, November 1986. 

G. D. Cohen, "A nonconstructive upper bound on covering radius", IEEE Trans. Inform. Theory, vol. IT- 
29, no. 3, pp. 352-353, May 1983. 

V. M.. Blinovskii, "A lower bound on the number of words of a linear code in an arbitrary sphere with 
given radius in F™" (in Russian), Probl. Pered. Inform. (Prob. Inf. Transm.), vol. 23, no. 2, pp. 50-53, 
1987. 



40 



[17] J. Chen, Da-Ke He, A. Jugmohan, "Achieving the rate-distortion bound with linear codes", IEEE Inform. 
Theory Workshop 2007, pp. 662-667, Lake Tahoe, California. 

[18] A. Barg and G. D. Forney Jr., "Random codes: Minimum distances and error exponents", IEEE Trans. 
Inform. Theory, vol. 48, no. 9, pp. 2568-2573, September 2002. 

[19] B. A. Nazer and M. Gastpar, "Computation over Multiple- Access Channels," IEEE Transactions on In- 
formation Theory, vol. 53, no. 10, pp. 3498-3516, Oct. 2007. 

[20] T. Philosof, A. Kishty, U. Erez and R. Zamir, "Lattice Strategies for the Dirty Multiple Access Channel" , 
Proceedings of IEEE International Symposium on Information Theory, July 2007, Nice, France. 

[21] D. Slepian, "Group codes for the Gaussian channel", Bell Syst. Tech. Journal, 1968. 

[22] G. D. Forney, Jr., "Geometrically uniform codes", IEEE Trans. Inform. Theory, vol. 37, no. 5, pp. 1241- 
1260, September 1991. 

[23] E. Biglicri and M. Elia, "On the existence of group codes for the Gaussian channel", IEEE Trans. Inform. 
Theory, vol. 18, no. 3, pp. 399-402, May 1972. 

[24] G. D. Forney, Jr. and M. D. Trott, "The dynamics of group codes: State spaces, Trellis diagrams, and 
Canonical encoders", IEEE Trans. Inform. Theory, vol. 39, no. 9, pp. 1491-1513, September 1993. 

[25] H. A. Loeliger, "Averaging bounds for lattices and linear codes", IEEE Trans. Inform. Theory, vol. 43, 
no. 6, pp. 1767-1773, November 1997. 

[26] H. A. Loeliger and T. Mittelholzer, "Convolutional codes over groups" , IEEE Trans. Inform. Theory,vo\. 42, 
no. 6, pp. 1660-1686, November 1996. 

[27] V. V. Vazirani, H. Saran and B. S. Rajan, "An efficient algorithm for constructing minimal trellises for 
codes over finite abelian groups", IEEE Trans. Inform. Theory, vol. 42, no. 6, pp. 1839-1854, November 
1996. 

[28] H. A. Loeliger, "Signal sets matched to groups" , IEEE Trans. Inform. Theory, vol. 37, no. 6, pp. 1675-1682, 
November 1991. 

[29] S. D. Berman, "On the theory of group codes", Kibernetika, vol. 3, no. 1, pp. 31-39, 1967. 

[30] G. D. Forney, Jr., "On the Hamming distance properties of group codes", IEEE Trans. Inform. Theory, 
vol. 38, no. 6, pp. 1797-1801, November 1992. 

[31] E. Biglicri and M. Elia, "Construction of linear block codes over groups", Proc. IEEE Interational Sympo- 
sium on Inform. Theory, San Antonio, TX, 1993. 

[32] J. C. Interlando, R. Palazzo and M. Elia, "Group block codes over nonabelian groups are asymptotically 
bad", IEEE Trans. Inform. Theory, vol. 42, no. 4, pp. '1277-1280, July 1996. 

[33] R. M. Tanner, D. Sridhara and T. Fuja, "A class of group structured LDPC codes", Proc. of ISCTA, 
Ambleside, 2001. 



41 



[34] G. Como and F. Fagnani, "The capacity of abelian group codes over symmetric channels", Submitted for 
publication. 

[35] F. Garin and F. Fagnani, "Analysis of serial turbo codes over abelian groups for geometrically uniform 
constellations" , Submitted for publication. 

[36] R. Ahlswcdc, "Group codes do not achieve Shannon's channel capacity for general discrete channels" , The 
Annals of Mathematical Statistics, vol. 42, no. 1, pp. 224-240, February 1971. 

[37] R. Ahlswede and J. Gemma, " Bounds on algebraic code capacities for noisy channels I", Information and 
Control, pp. 124-145, 1971. 

[38] R. Ahlswede and J. Gemma, "Bounds on algebraic code capacities for noisy channels II" , Information and 
Control, pp. 146-158, 1971. 

[39] J. Korner and K. Marton, "How to encode the modulo-two sum of binary sources," IEEE Trans. Inform. 
Theory, vol. IT-25, pp. 219-221, March 1979. 

[40] A. D. Wyner, "Recent results in Shannon theory," IEEE Trans, on Inform. Theory, vol. 20, pp. 2-10, 
January 1974. 

[41] A. D. Wyner, "On source coding with side information at the decoder," IEEE Trans, on Inform. Theory, 
vol. IT-21, pp. 294-300, May 1975. 

[42] R. Ahlswede and J. Korner, "Source coding with side information and a converse for degraded broadcast 
channels," IEEE Trans. Inform. Theory, vol. IT- 21, pp. 629-637, November 1975. 

[43] H. Viswanathan, Z. Zhang, and T. Berger, "The CEO problem," IEEE Trans, on Inform. Theory, vol. 42, 
pp. 887-902, May 1996. 

[44] A. D. Wyner and J. Ziv, "The rate-distortion function for source coding with side information at the 
decoder," IEEE Trans. Inform. Theory, vol. IT- 22, pp. 1-10, January 1976. 

[45] T. Berger, "Multiterminal source coding," in Lectures presented at CISM summer school on the Inform. 
Theory approach to communications, July 1977. 

[46] S.-Y. Tung, Multiterminal source coding. PhD thesis, School of Electrical Engineering, Cornell University, 
Ithaca, NY, May 1978. 

[47] T. Berger and R. W. Yeung, "Multiterminal source coding with one distortion criterion", IEEE Trans, on 
Inform. Theory, vol. IT-35, pp. 228-236, March 1989. 

[48] R. Dobrushin and B. Tsybakov, "Information transmission with additional noise," IRE Trans. Inform. 
Theory, vol. IT- 18, pp. S293-S304, 1962. 

[49] H. S. Witsenhausen, "Indirect rate distortion problems," IEEE Trans. Inform. Theory, vol. IT- 26, pp. 518- 
521, September 1980. 



42 



[50] S. Gelfand and M. Pinsker, "Coding of sources on the basis of observations with incomplete information," 
Problemy Peredachi Informatsii, vol. 15, pp. 45-57, Apr- June 1979. 

[51] T. Berger, "Multiterminal source coding", Lecture notes presented at the 1977 CIS M summer school, Udine, 
Italy, pp. 171-231, July, 1977. 

[52] R. Ahlswede and T. S. Han, "On source coding with side information via a multiple-access channel and 
related problems in multi-user information theory," IEEE Trans, on Inform. Theory, vol. 29, pp. 396-412, 
May 1983. 

[53] I. Csiszar, "Linear codes for sources and source networks: Error exponents, universal coding," IEEE Trans. 
Inform. Theory, vol. IT- 28, pp. 585-592, July 1982. 

[54] D. Krithivasan and S. S. Pradhan, "Lattices for distributed source coding: Jointly Gaussian Sources and 
Reconstruction of a linear function," Submitted to IEEE Trans. Inform. Theory. 

[55] R. Zamir, S. Shamai, and U. Erez, "Nested linear/lattice codes for structured multiterminal binning," 
IEEE Trans. Inform. Theory, vol. LT-48, pp. 1250-1276, June 2002. 

[56] S. Shamai, S. Verdu, and R. Zamir, "Systematic lossy source/channel coding," IEEE Trans, on Inform. 
Theory, vol. 44, pp. 564-579, March 1998. 

[57] I. Csiszar and J. Korner, Information Theory: Coding Theorems for Discrete Memoryless Systems. Aca- 
demic Press Inc. Ltd., 1981. 

[58] T. S. Han and K. Kobayashi, "A dichotomy of functions F(X,Y) of correlated sources (X,Y)," IEEE Trans, 
on Inform. Theory, vol. 33, pp. 69-76, January 1987. 

[59] D. S. Dummit and R. M. Footc, Abstract Algebra. John Wiley & sons Inc., 2004. 

[60] A. G. Kurosh, The Theory of Groups. Chelsea publishing company, 1960. 

[61] S. Janson, "New versions of Sucn's correlation inequality," Random Structures Algorithms, vol. 13, pp, 467- 
483, 1998. 

[62] D. Slepian and J. K. Wolf, "Noiseless coding of correlated information sources, " IEEE Trans, on Inform. 
Theory, vol. 19, pp. 471-480, July 1973. 

[63] A. D. Wyner, "On source coding with side information at the decoder," IEEE Trans, on Inform. Theory, 
vol. 21, pp. 294-300, May 1975. 

[64] W. Gu, S. Jana and M. Effros, "On approximating the rate regions for lossy source coding with coded and 
uncoded side information", Proc. IEEE International Symposium on Inform. Theory, Toronto, Canada, 
2008. 

[65] S. Boyd and L. Vandenberghe, Convex Optimization. Cambridge University Press, 2004. 



43 



[66] M. V. Eyuboglu and G. D. Forney, "Lattice and trellis quantization with lattice and trellis-bounded 
codebooks-High rate theory for memoryless sources," IEEE Trans. Inform. Theory, vol. 39, pp. 46-59, 
Jan. 1993. 

[67] T. Philosof and R. Zamir, "The rate loss of single-letter characterization: The "Dirty" multiple access 
channel," Available at \arXiv:0803.1120b 3 



44 



